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1.0  INTRODUCTION 


Solid  rocket  propellants,  being  highly  energetic  compounds,  may, 
upon  exposure  to  suitable  stimulus,  deflagrate  or  detonate.  For  rocket 
propulsion,  these  energetic  compounds  must  deflagrate  in  a predictable 
and  controllable  manner,  converting  the  stored  potential  chemical  energy 
into  high  temperature  gases  which,  when  exhausted,  result  in  a thrust 
imparted  to  the  missile.  Unfortunately,  this  same  stored  chemical 
energy  can  also  be  released  as  a detonation.  Usually,  the  stimulus 
causing  detonation  is  provided  by  mechanical  action  and  is  referred  to 
as  shock-to-detonation  initiation  (SDI)  (analogous  to  initiation  by  a 
small  amount  of  primary  explosive).  However,  detonation  can  also  occur 
as  a transition  from  deflagration  without  any  external  mechanical  initi- 
ating source.  In  this  case,  called  deflagration-to-detonation  transition 
(DDT),  the  initiating  energy  originates  in  the  deflagration  itself  as  a 
result  of  the  synergistic  interaction  between  pressure  and  mass  burning 
rate,  leading  to  abnormally  high  rates  of  pressure  rise.  This  transition 
is  thought  to  occur  as  follows: 


The  first  transition  might  be  thought  of  as  the  change  from  normal 
conductive  (layerwise)  combustion  to  the  more  rapid  abnormal  convective 
combustion  (burning  "within"  the  sample).  Pressure  generation  rate  is  a 
function  of  the  mass  burning  rate,  and  mass  burning  rate  is  proportional 
to  the  product  of  linear  burning  rate  (surface  regression  rate)  and  the 
burning  surface  area.  Since  the  known  dependence  of  surface  regression 
rates  on  pressure  cannot  alone  account  for  the  pressurization  rates 
required  to  cause  transition  to  detonation,  then  increased  burning 
surface  area  must  also  be  involved  in  producing  the  required  mass  burning 
rates,  hence  pressurization  rates.  Two  requirements  must  be  satisfied 
in  order  to  transit  from  normal  conductive  burning  to  the  abnormal 
convective  burning:  (1)  the  existence  of  extra  surface  area  (defects) 
prior  to  combustion  or  its  creation  during  combustion,  and  (2)  the 
development  of  combustion  on  the  increased  surfaces.  The  correctness  of 
the  above  seems  to  be  confirmed  by  the  lack  of  experimental  evidence 
showing  a consolidated  charge  undergoing  DDT.  In  fact,  investigators 
studying  DDT  phenomena  often  "shred"  their  samples. 
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Given  the  simple  flow  chart  presented  above  and  the  two  criteria 
outlined  in  the  preceding  paragraph,  then  one  of  the  most  important 
areas  to  be  considered  in  trying  to  understand  DDT  is  the  penetration 
and  establishment  of  the  combustion  front  in  propellant  defects.  It  is 
the  purpose  of  this  survey  to  present  what  is  known  of  the  early  develop- 
ment of  combustion  into  defects  - a necessary  condition  for  the  later 
transition  to  full  convective  burning  and  detonation. 

Twc  aspects  of  combustion  are  considered:  the  onset  of  anomalous 
combustion  and  the  progression  of  combustion  into  the  defects.  The 
first  is  primarily  qualitative  and  relates  to  experimental  observation 
of  a singularity  indicating  whether  or  not  convective  burning  occurred 
(go-no  go  type  of  test).  Results  of  onset  experiments  are  commonly 
displayed  as  a graphical  relation  which  shows  combinations  of  pressure 
and  defect  dimension  for  which  combustion  either  does  or  does  not  pene- 
trate the  defect.  The  coordinates  of  such  plots  are  usually  referred  to 
as  critical  pressure  or  critical  defect  size  (crack  width  or  pore  diam- 
eter); it  is  actually  the  relation  between  the  two  which  is  critical, 
and  no  importance  should  be  attributed  to  the  variable  which  is  so 
characterized.  Other  terms  synonymous  with  critical  pressure  are: 
breakdown  pressure,  transition  pressure,  threshold  pressure,  and  pressure 
at  which  stable  burning  is  impaired.  The  critical  pore  size  is  sometimes 
referred  to  as  the  threshold  crack  width  or  pore  diameter. 

The  second  aspect  of  combustion  is  more  quantitative,  is  not  as 
extensively  investigated,  requires  more  sophisticated  measurement  tech- 
niques, and  yields  information  on  the  rate  of  propagation  of  the  ignition 
front  into  the  propellant  defect  up  to,  but  not  including,  generation  of 
weak  shocks. 

At  this  point  it  is  well  to  divert  the  discussion  to  define  the 
terminology  relevant  to  the  various  rates  (velocities)  of  processes 
encountered  in  the  survey. 

Conductive  rate.  This  is  the  rate  of  regression  of  a solid  surface 
of  propellant  in  chi/s.  It  is  often  called  "normal”  burning  rate,  surface 
regression  rate,  o^r  simply  burning  rate.  The  latter  term  will  not  be 
used  in  this  survey  , because  of  possible  confusion  with  mass  burning 
rate.  ' 

Ignition  rate.  This  is  the  rate  in  cm/s  at  which  a state  of  igni- 
tion (self  sustained  combustion)  is  propagated  into  a single  pore  or 
porous  bed  of  propellant.  Synonymous  expressions  are  ignition  propaga- 
tion rate,  convective  burning  rate,  convective  propagation  rate,  or 
convective  rate. 
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Flow  rate.  This  is  the  rate  in  cm/s  at  which  hot  gases  flow  in  a 
single  pore  or  porous  bed,  with  or  without  ensuing  ignition.  As  will  be 
described  in  the  appropriate  sections  of  the  survey,  the  gases  always 
flow  ahead  of  the  ignition  front  during  early  stages  of  development  of 
convective  burning. 

Mass  burning  rate.  This  is  the  rate  of  generation  of  combustion 
products  ill  g/cm^-s.  As  stated  earlier,  it  is  determined  both  by  the 
conductive  rate  and  the  burning  surface  area,  and  is  the  quantity  direct- 
ly responsible  for  rate  of  pressure  rise. 

Much  of  the  experimental  and  theoretical  work  on  the  early  stages 
of  convective  burning  has  been  conducted  over  the  last  30-A0  years  in 
the  Soviet  Union.  The  work  covers  a broad  spectrum  of  experimental 
conditions,  sample  geometries,  and  propellant  types.  It  is  important  to 
note  that  the  term  propellant,  when  used  to  describe  Soviet  work,  does 
not  imply  an  operational  rocket  propellant,  but  is  used  synonymously 
with  such  expressions  as  propellant  powder,  explosive,  porous  system, 
porous  powder,  porous  charge,  and  others.  Regardless  of  the  name,  a 
common  feature  of  the  samples  used  by  the  Soviets  in  porous  bed  experi- 
ments is  that  they  are  prepared  by  pressing  or  pouring  mixtures  of 
ingredients  in  powdered  form.  While  the  ingredients  may  include  poly- 
meric materials,  no  use  is  madr  of  the  binding  properties  of  the  polymer 
to  prepare  a monolithic  sample. 

This  report  is  primarily  a survey  of  the  Soviet  literature  but  also 
includes  non-Soviet  literature  relevant  to  early  convective  burning. 

The  purpose  is  to  bring  together  in  an  organized  fashion  the  available 
information  on  the  subject  to  serve  as  a point  of  departure  for  future 
experimental  and  theoretical  work.  In  addition,  answers  were  being 
sought  to  several  broad  questions. 

1.  What  propellant  characteristics  are  responsible  for  the  onset 
and  development  of  the  early  steps  of  convective  burning? 

2.  To  which  characteristics  are  this  onset  and  development  most 
and  least  sensitive? 

3.  What  tests  and  test  conditions  are  suitable  for  determining 
these  characteristics? 

4.  What  tests  and  test  results  ir.ost  nearly  simulate  operational 
firings  of  rocket  motors? 

The  next  section  of  the  report  explains  the  organization  o*'  the 
survey . 
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2.0  SURVEY  ORGANIZATION 


The  literature  describing  the  early  stages  of  convective  burning 
consists  of  both  experimental  (Section  3)  and  theoretical  (Section  4) 
approaches  to  the  problem.  Two  separate  sections  were  chosen  because, 
unfortunately,  there  is  little  systematic  interaction  between  the  two 
approaches . 


2.1  ORGANIZATION  OF  SECTION  3:  RESULTS  FROM  EXPERIMENTS 

There  are  many  aspects  of  convective  burning,  and  experiments  have 
been  designed  to  understand  various  combinations  of  variables.  The  many 
variables  to  be  controlled  can  be  grouped  as  follows. 

2.1.1  Sample 

Composition  (homogeneous  and  heterogeneous)  and  sample  geometry 
have  been  widely  studied.  Two  main  geometric  divisions  are  single  pores 
(or  channels  bounded  by  propellant)  and  porous  beds  of  propellant.  A 
third  type  of  sample,  damaged  propellant,  was  studied,  but  is  not  in- 
cluded in  the  morphological  organization  because  of  the  limited  amount 
of  work  conducted.  The  first  is  characterized  by  known  length  and 
cross-sectional  shape.  The  second  is  an  array  of  pores  which  are  numer- 
ous enough  that  they  must  be  described  by  statistical  average  quantities 
such  as  porosity,  pore  size  distribution,  particle  size,  specific  surface 
area,  and  permeability  to  fluid  flow.  Methods  of  measuring  these  quan- 
tities as  well  as  theories  of  flows  through  packed  beds  are  presented  in 
numerous  standard  references  [Scheidegger  (1974),  Muskat  and  Wyckoff 
(1946),  Bear  (1972),  Carman  (1956)]  as  well  as  in  one  of  the  surveyed 
reports  [Belyaev  (1973),  Sections  4 and  5].  Two  comments  on  porous 
materials  are  appropriate  at  this  point.  First,  the  Russian  theoretical 
approaches  to  early  convective  burning  are  based  entirely  on  the  perme- 
ability concept  rather  than  on  the  fluidized  bed  concept  used  in  many 
western  models  of  the  development  of  shock  waves.  Second,  the  single 
pore  and  the  statistical  porous  bed  represent  two  extremes  with  the 
actual  conditions  of  irregularly  branched  cracks  in  between  where  no 
precise  quantification  of  the  defect  has  yet  been  devised  for  correlation 
with  experimental  results. 

Both  of  the  m.ain  types  (single  pore  vs.  porous  bed)  may  be  further 
characterized  to  include: 

1.  Tile  iMubeddod  or  built-in  charge.  In  Liiis  case  the  sample  is 
eoiii])  1 e I e'l  y eiicaia'd  in  pas  impi'rmeable  material  (Figure  la,  b)  assuring 
that  the  defect  is  not  prepressurized  by  the  gases  that  would  pressurize 
a closed  vessel.  This  condition  simulates  rocket  cliarge  defects  that 
are  not  orig.inally  opiai  to  the  bore,  i.c.,  the  burning  surface  must 
regj-i.ss  to  the  defec't. 
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FIGURE  1.  Test  Sample  Geometries. 


2.  The  open  defect.  Here  the  defect  is  left  open  to  the  bomb 
pressure  and  the  ignition  event  at  the  ignited  end.  This  case  has 
several  subcases: 


a.  The  unignited  end  is  sealed  (the  so-called  blind  or 
closed  defect  shown  in  Figure  Ic,  d) . 

b.  The  unignited  end  is  open  to  the  bomb  pressure  (Figure 
le,  f) . 

c.  The  unignited  end  is  connected  through  a valve  to  a 
pressure  less  than  the  bomb  pressure  (1  atm  [0.1  MPa]  was 
used  by  Krasnov)  (Figure  Ig) . 

2.1.2  Environment  and  Stimulus 


In  general,  convective  burning  is  studied  experimentally  by  burning 
one  of  the  propellant  samples  discussed  above  under  controlled  condi- 
tions. Development  of  convective  burning  depends  upon  a pressure  differ 
ential  between  the  external  pressure  and  the  pressure  within  the  defect. 
The  external  pressure  may  be  provided  by  carrying  out  the  experiment 

(1)  in  a closed  bomb  with  only  a slight  pressure  rise  during  the  test, 

(2)  in  a closed  bomb  in  which  the  pressure  rises  rapidly  and  significant 
ly  during  the  test,  or  (3)  at  ambient  atmospheric  pressure.  In  the 
literature  reviewed,  the  rapid  pressure  rise  is  provided  (with  a single 
exception)  by  combustion  of  the  test  charge  in  a high  loading  density 
combustion  bomb  (HLDCB)  with  a typical  loading  density  (grams  of  propel- 
lant per  cubic  centimeter  of  bomb  volume)  of  0.02-0.1  g/cm^.  Tests 
leading  to  small  pressure  rises  are  conducted  in  a low  loading  density 
combustion  bomb  (LLDCB)  with  a typical  loading  density  less  than  0.0002 
)’/cm-^,  where  the  contribution  of  gases  from  the  burning  charge  are 
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negligible  compared  with  the  initial  gas  pressure.  The  HLDCB  (Russian 
equivalent:  manometer  bomb  or  bomb  of  rising  pressure)  and  LLDCB 
(Russian  equivalent:  bomb  of  constant  pressure)  were  discussed  by  Boggs 
(1976).  In  the  single  exception,  noted  above,  the  high  rate  of  pressure 
rise  is  provided  by  admitting  nitrogen  gas-  through  a regulator  to  the 
bomb  during  the  test.  It  should  be  noted  that  the  quantitative  effect 
of  rate  of  pressure  rise  in  the  HLDCB  has  not  been  studied  extensively. 
The  pressure  in  the  LLDCB  also  rises  during  a test  (the  magnitude  and 
rate  may  be  partially  offset  by  using  a surge  tank);  the  few  atmospheric 
tests  reported  eliminate  this  effect.  There  are  no  tests  employing 
pressures  above  atmospheric  with  strictly  zero  pressure  rise  rate.  Such 
capability  would  be  useful  in  exploring  dynamic  pressure  effects  but 
would  not  simulate  any  practical  condition. 

While  the  pressure  in  the  HLDCB  is  provided  primarily  by  sample 
combustion,  the  pressure  in  the  LLDCB  must  be  adjusted  to  the  desired 
operating  point  by  introducing  a gas,  usually  nitrogen  or  argon,  from  an 
external  source.  The  implications  of  this  difference  in  operating  mode 
are  considerable.  In  addition  to  requiring  many  more  tests  to  span  a 
given  pressure  range  [Boggs  (1976)]  the  driving  force  causing  penetration 
of  combustion  into  propellant  defects  may  be  quite  different  in  the  two 
cases.  In  the  HLDCB,  the  forcing  pressure  differential  arises  because, 
at  the  beginning  of  the  test,  the  pressure  external  to  the  defect  rises 
faster  than  the  internal  pressure.  In  the  LLDCB,  there  are  two  possibil- 
ities. First,  if  the  test  sample  is  "embedded,”  so  that  its  pores  (at 
low  pressure,  e.g.,  1 atm  [0.1  MPa])  are  sealed  against  the  initial  bomb 
pressure  by  an  impermeable  layer  of  auxilliary  propellant,  then  at  the 
instant  of  burnthrough  of  this  layer,  the  pressure  differential  is  set 
by  the  difference  between  pore  pressure  and  the  bomb  pressure.  Obviously 
this  arrangement  is  feasible  only  for  the  blind  pore  or  porous  bed. 
Second,  if  the  test  sample  is  open  to  the  bomb  pressure,  the  pores  are 
initially  filled  with  cold,  high  pressure  inert  gas  used  to  pressurize 
the  bomb.  Now  the  pressure  differential  responsible  for  convective 
burning  originates  in  the  slow  pressure  rise  in  the  LLDCB  or  in  the 
combustion  process  itself  (dynamic  pressure),  which  generates  a slightly 
higher  pressure  at  the  surface  of  the  propellant  than  in  the  surrounding 
volume.  Moreover,  the  presence  of  the  cold  gases  in  the  pore  makes 
penetration  of  combustion  into  the  pores  more  difficult  because  the 
initial  gases  must  be  displaced  or  heated  before  ignition  of  the  pores 
can  occur. 

The  morpliology  described  above  is  graphically  presented  in  Table  1, 
which  also  indicates  references  where  relevant  information  was  obtained 
for  this  survey.  The  columns  of  this  matrix  are  concerned  with  the 
sample.  This  broad  area  is  divided  into  single  pore  and  porous  beds. 

The  single  pore  is  further  divided  into  embedded  (Figure  la),  blind 
(Figure  Ic) , macaroni  (open-open  as  shown  in  Figure  le),  and  valved 
(Figure  1 g)  . The  ijorous  bed  is  divided  into  embedded  (Figure  lb),  blind 
(Fig.urt:  Id),  and  open-open  (Figure  If).  The  rows  in  the  matrix  are 


8 


TABLE  1.  Morphology  of  Literature  on  Experimental  Convective  Burning 


[ 

i 

} 


i 


N\;C  TP  6007 


.i 

1 

i 


9 


NWC  TP  6007 


concerned  with  the  environment  and  stimulus  provided  by  the  test  ap- 
paratus. The  rows  include;  (1)  the  p > 1 atm  (0.1  MPa),  p w 0 condi- 
tions associated  with  LLDClis  such  as  a Crawford  bomb,  (2)  the  p > 1 atm 
(0.1  MPa) , p>>0  conditions  associated  with  HLDCBs,  and  (3)  a series  of 
tests  where  the  environment  was  1 atm  (0.1  MPa)  and  no  imposed  pressure 
change  (stimulus)  was  applied. 

The  matrix  is  further  divided  by  whether  the  onset  (0)  or  develop- 
ment (D)  of  convective  combustion  was  studied. 

Section  3 (Experimental  Studies)  will  follow  the  organization  of 
Table  1.  That  is,  each  column  will  be  discussed  in  order  in  its  entire- 
ty. 


2.2  ORGANIZATION  OF  SECTION  4:  RESULTS  FROM  THEORETICAL  APPROACHES 

There  are  two  kinds  of  theoretical  models  dealing  with  DDT  pro- 
cesses. The  first,  more  complicated  approach,  is  based  on  conservation 
laws,  expressed  as  differential  equations,  augmented  by  empirical  rela- 
tionships for  friction,  heat  transfer,  and  propellant  ignition  and 
combustion  [Kuo  (1977);  Takata  and  Wiedermann  (1976);  Kuo  (1973); 

Krier  and  Gokhale  (1976)].  An  advantage  is  that,  assuming  a correct 
analysis,  given  proper  physical  data,  and  using  a high-speed  computer, 
one  may  calculate  meaningful  detailed  solutions  representing  growth  and 
propagation  of  combustion  in  single  channels  or  in  porous  beds.  There 
are  two  disadvantages  in  such  an  attack.  First,  solutions  do  not  exhibit 
transition  from  conductive  to  convective  burning;  hence  no  onset  criter- 
ion can  be  established  from  solutions.  Second,  computations  are  expen- 
sive, so  that  it  is  not  economical  to  carry  out  parameter  and  sensitivity 
studies.  As  a result,  detailed  studies  of  trends  of  solutions  with 
relevant  changes  in  input  data  have  not  been  made. 

By  contrast,  the  other,  simpler  approach  is  characterized  by 
several  approximations  in  order  to  reduce  complexity  and  permit  solutions 
to  be  obtained  without  the  need  for  extensive  computer  usage.  While  the 
results  are  not  quantitatively  accurate,  they  usually  reveal  qualitative 
trends  which  may  be  compared  to  experimental  results  if  it  is  realized 
that  detailed  agreement  is  not  to  be  expected.  It  is  the  literature 
concerning  these  simplified  modeling  efforts  which  is  reviewed  in  Section 
4 of  this  report. 

It  is  convenient  to  separate  the  theoretical  literature  into  major 
sections  treating  single  pores  and  porous  beds,  as  was  done  in  the 
experimental  section.  Here  the  similarity  ends  for  throe  reasons: 

(1)  the  scope  and  extent  of  theoretical  literature  is  far  more  limited 
so  that  further  categorization  similar  to  that  outlined  in  Section  2.1 
is  not  warr/mted;  (2)  emphasis,  in  the  case  of  simi)lificd  theories  (as 
opposed  to  eomplex),  i ,s  oriented  more  toward  stages  in  the  deve  lopiiu'n  t 
et  convective  burning  tiian  toward  satiiivLe  houndarv  conditions  and  pressure 
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environment;  and  (3)  there  are  miscellaneous  theoretical  analyses  which 
do  not  even  fit  the  experimental  morpliology.  The  emphasis  on  stages  is 
really  part  of  the  simplified  approacli  in  which  tEie  overlapping  (in 
time)  pEienomena  are  heuristically  decoupled  in  order  to  enable  a tract- 
able solution  of  tlie  problem. 


Tliere  are  three  requirements  for  the  breakdown  of  stable  conductive 
(normal)  burning  in  a single  pore  or  porous  bed.  First,  the  hot  combus- 
tion gases  must  flow  into  tEic  defect  at  a velocity  greater  than  the 
conductive  regression  rate.  Otherwise  the  defect  would  not  exist  for 
the  gases  to  penetrate.  Second,  the  hot  gases  must  heat  the  pore  walls 
to  a condition  of  self-sustained  combustion . Third,  this  condition 
(usually  called  ignition)  must  continue  to  propagate  into  tEie  single 
pore  or  porous  bed  at  a velocity  exceeding  the  conductive  regression 
rate.  The  first  requirement  poses  a hydrodynamic  problem  for  whicli 
simplified  solutions  liave  been  found  for  single  pores  and  porous  beds. 

In  all  instances  the  energy  equation  is  either  ignored  or  approximated 
by  an  assumption  of  isothermal  flow.  Ttie  second  requirement  Involves 
heat  transfer  and  some  kind  of  ignition  criterion,  both  of  which  are 
included  in  a non-mechanistic  manner.  In  the  literature  surveyed,  only 
single  pores  have  been  dealt  with.  The  requirement  for  propagation  of 
ignition  is  presented  in  only  one  paper,  again  limited  to  a single  pore. 
The  extension  of  single  pore  analyses  to  porous  bed  situations  and  vice 
versa  is  not  exploited  to  any  great  degree  because  of  the  difficulty  and 
uncertainty  of  assigning  an  effective  pore  size  to  a porous  bed.  The 
attempt  was  made,  in  connection  with  the  concept  of  Andreev  number,  by 
Margolin  and  Chuiko  (1966)  and  by  Bobolev  (1966).  In  botli  references, 
only  the  conditions  for  onset  of  convective  burning  were  considered; 
extension  of  single  pore  analyses  of  ignition  and  propagation  of  ignition 
to  porous  beds  is  not  made. 


Section  4 of  this  report  presents  a review  of  the  available  theo- 
retical literature  which  treats  the  three  requirements  for  convective 
burning  in  single  pores  and  porous  beds.  In  addition,  mention  is  made 
of  several  miscellaneous  papers,  dealing  with  melt  layers,  oscillatory 
pressure,  and  dynamic  (combustion  generated)  pressure. 


3.0  RieSULTS  I- ROM  EXPERIMENTS 


The  outline  and  organization  of  tills  section  were  presented  in 
Section  2.1. 

3.1  EXPERIMENTS  USING  THE  SINGE, E PORE  (OR  SINGLE  CHANNI'L)  GEOMETRY 

3.1.1  Embedded  Single  refect  Tested  in  the  LLDCB 

l{ccalling  the  discussion  of  using  the  I.LDCB  to  study  combustion  of 
samples  having  single  defects  (Section  2.1),  two  possibilities  were 
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mentioned.  If  the  test  sample  was  embedded  (Figure  la),  the  defect 
(presumably  having  1 atm  (0.1  MPa)  pressure)  was  sealed  against  the 
initial  bomb  pressure  by  a layer  of  auxilliary  propellant.  At  the 
instant  of  burnthrough  of  this  layer,  the  defect  would  "see"  the  much 
higher  pressure  of  the  bomb,  and  the  required  pressure  differential 
would  be  established. 

Only  limited  experimental  work  has  been  carried  out  using  embedded, 
single  defects  in  the  LLDCB  (Godai  (1970)].  The  defects  consisted  of 
narrow  slits,  25  mm  in  length,  constructed  by  assembling  two  5 X 5 mm 
strands  of  propellant  with  an  intervening  gap  varying  from  0.05  to  0.5 
mm.  The  propellant  was  a formulation  consisting  of  20Z  polyester  and 
80%  trimodal  AP  (nominal  particle  sizes:  24,  160,  and  480  urn;  no  propor 
tions  given).  Experiments  were  conducted  in  nitrogen  and  air  at  pres- 
sures up  to  4 atm  (0.4  MPa)  with  no  noticeable  effect  of  the  gas.  It 
was  stated  that  tests  run  in  helium  gave  sliglitly  different  results. 

Only  the  onset  conditions  were  obtained  as  shown  in  Figure  2.  Although 
the  pressure  range  is  quite  limited,  results  are  consistent  with  the 
trend  chat  threshold  pressure  increases  as  defect  size  decreases. 


I’RFSSURf.  p.  MPa 

FIGURE  2.  Burning  Rate  and  Threshold 
Crack  Width  for  l.ow  Pressun-  Convective 
Burning  of  PotyesLer/AP  Propellant. 

(Data  taken  from  Godai  [1970],  Figure  5.) 


NWC  TP  6007 


3 . 1 . 2 . 1 The  Blind  Slnj^le  Detect  in  tiu'  LLUCP 

In  this  configuration.  Figure  Lc,  the  def(ut  of  the  test  sanipli-  is 
open  to  tlie  bomb  atmosphere  and  is  intially  filled  with  tlie  cold,  high 
pressure  inert  gas  used  to  pressurize  the  bomb.  The  presence  of  the 
cold,  inert  gases  in  the  defects  makes  penetration  of  combustion  into 
the  defects  more  dtlfieult  bi-cause  the  initial  gases  must  be  displaced 
or  heated  before  ignition  can  occur.  This  is  espeeialJy  true  for  the 
blind  defect  because  the  only  exit  for  gas  flow  (the  mouth  of  the  defect) 
Is  the  entrance  that  tlie  combustion  front  must  penetrate.  T!ie  slow 
pressure  rise  associatt'd  witli  tlie  LLDt^B  compounds  the  diflicuity  ol 
combustion  penetrating  tin?  defect  sinci'  the  pressure  differential  respon- 
sible for  causing  convective  burning  comes  entirely  from  tlte  combu.sLion 
process  itself. 


Results  presented  by  Bely.aev  (1973)  (Section  22)  piovide  limited 
information  on  tlie  rate  of  spontaneous  propagation  of  combustion  into 
blind  flat  cracks  for  conditions  in  whch  only  a small  pre.ssure  ri.sc 
occurs  in  the  crack.  Two  typiis  of  cracks  were  used:  (1)  two  opposing 
"composite  propellant"  (no  composition  or  roniiuLatioils  given)  surlaces 
and  (2)  opposing  "propellant"  and  plexiglass  surfacivs.  The  qualitative 
effect  of  pre.ssure  and  crack  width  on  convective  rates  is  sliown  in 
Figure  3.  The  ordinate  scale  is  tlie  ratio  of  convective  to  conductive 
rate  at  the  ' ressure  Involved.  It  Is  seen  that  the  propagation  rate 
Increases  wi, h increasing  pressure  and  with  decreasing  crack  width 
until,  for  a minimum  crack  widtii,  propagation  no  longer  occurs.  Detailed 
studies  of  the  region  to  the  right  o!  the  maxima  of  Figure  3 are  shown 
in  Figure  4 for  a composite  propellant.  It  is  not  clear  whieh  crack 
type  was  used  in  obtaining  Figure  4 anti  Eq . (1).  Tlie  data  arc  repre- 

sented by  the  expression: 

- / t “0-  ^ n >1  I “0.  h ( , t 

(v . /r)w  p =0.26  cm  ;itm  (J1 

ig 

where 


ig 

r 

w 

p 


ignition  propagation  rate 

conductive  rate 
crack  widtli 
pressure 


Although  it  is  claimed  that  a similar  result  cun  be.  derived  t rom  theory, 
no  such  derivation  or  reference  to  it  was  given. 


Without  actual  supiiortive  data,  the  following  conclusions  are 
stated  (liely.iev  [1973,  Section  22])  regarding  ])ropagation  rates  into 
blind  channels  at  const  .iiit  jiressure: 


13 


NUC  TP  6007 


fLAT  CRACK 


FIGURE  3.  Quaiitat.iviJ  Effect  of  Crack.  Width 
and  Pressure  on  Convective  Burning  Rate. 
[Data  taken  from  Belyaev  (1973),  Figure  51.] 


FIGURE  4.  Effect  of  Crack  Width  and  Pres.sure 
on  Convective  Burning  Rate  of  Composite  Pro- 
pellant Powder.  (Data  taken  from  Belyaev 
(1973) , Figure  52.  ] 
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1.  Ij'nition  rate  increases  wltli  channel  length  at  constant 
length- to-width  ratio. 

2.  Ignition  rate  is  greater  for  two  opposing  burning  .surfaces 
than  for  one  inert  and  one  burning 

3.  Ignition  rate  is  increased  by  channel  wall  roughness. 

4.  Ignition  rate  is  greater  in  a flat  crack  than  in  a circular 
channel  witli  equal  hydraulic  diameter. 

Others  have*  also  investigated  the  blind  single  pores  in  the  LLDCB 
but  this  was  only  one  of  the  configurations  and  apparatus  combinations 
studied  in  their  experiments.  Ratlier  tliau  discuss  the  results  here  - 
they  are  more  meaningful  in  context  - only  mention  of  the  references 
will  be  made  at  this  point.  Hobolev  (1965a)  i-ompared  threshold  condi- 
tions for  open  and  closed  flat  cracks  of  two  different  lengths  under 
conditions  of  constant  (LLDCB)  and  rising  (HLDCB)  pressure  (discussed  in 
Section  3. 1.3.2).  Margolin  and  Margulis  (1969)  studied  nitroglycerin 
powder  in  several  configurations  (discussed  in  Section  3.1.3.  1). 

Belyaev  (1969)  conducted  ti’sts  on  se(H)iidary  explosives,  nitroglycerin 
powder  and  AP  mixtures  in  flat,  closed  and  o)>en  cracks  of  varying  widths 
(discussed  in  Section  3. 1.3.1).  Prentice  (1962  and  1977)  lias  studied 
both  blind  and  open  single  pores  in  the  LLDCB  (discussed  in  Section 
3. 1.3.1).  Discussion  of  the  results  is  deferred  to  the  sections  indi- 
cated. 

3 . 1 . 2 . 2 Blind  Single  Defect  in  the  HLDCB 

The  HLDCB  was  used  to  determine  threshold  pressure  vs.  pore  width 
for  several  "ptopeilnnts"  in  a narrow  slit  configuration  [Belyaev 
(1973)].  Again,  two  types  of  slits  were  used.  As  described  in  Section 
3.1.2. 1,  results,  shown  in  Figure  5,  were  independent  of  the  configura- 
tion. This  is  no  contradiction  of  conclusion  2 of  the  preceding  section 
which  applies  to  LLDCB  and  ignition  rate.  As  with  many  of  tlu-  HLDCB 
tests  reported,  the  rate  of  change  of  pressure  was  not  recorded  nor  was 
its  effect  determined.  It  can  be  .scon  by  comparing  the  carves  presented 
that  of  the  two  composite  propellant  powders  studied,  the  one  with  the 
higher  burning  rate  develojis  convecLivt'  burning  at  a lower  pressure  for 
a given  crack  width.  Comparing  composite  and  NG  powders,  one  notes  that 
the  former  are  less  stable.  A possible  e.xplanation  is  the  smaller  flame 
standoff  distance  and  heterogeneity  in  the  flame  xone  owing  to  diffusion 
flames  of  the  composites,  which  al]f)ws  liotter  gases,  with  greater  ca])ac- 
ity  for  igniting  the  walls,  to  enter  the  crack.  Siir.]-)!!  fi  cation  of  an 
already  simple  theoretical  apprtiach  (see  Theoretical  section)  provides 
an  express!  in  connecting  pore  width  and  threshold  pressure; 

l+2n  2 


where 


= threshold  i^ressure 

n = pr('.ssure  exponent  in  conductive  rate  law 
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1 - RDX 

2 - Nitroglycerin  Powder 

3 - Slow  Burning  Composite  Propellant  Powder 

4 - Fast  Burning  Composite  Propellant  Powder 


FIGURE  5.  Dependence  of  ThreshoJ.d  Pressure  on 
Crack  Width  for  a Flat  Crack.  [Data  taken  from 
Belyaev  (1973),  Figure  38.) 


The  value  of  n is  unity  for  RDX  and  0.5  for  the  other  propellants;  good 
agreement  between  Eq . (2)  and  Figure  5 is  found  for  small  pore  widths. 

Except  for  RDX,  there  is  a break  in  tlie  vs.  w relationshi ' at  large 
crack  widths;  this  is  explained  as  resulting  from  one  of  the.  simp  .ifying 
assumptions  (constant  Nusselt  number)  made  In  deriving  Eq . (2).  Tliere 

is  no  explanation  offered  regarding  the  absence  of  this  feature  in  the 
data  for  RDX. 

The  work  of  Belyaev  (1973)  (Section  22)  on  blind  cracks  in  LLDCB 
was  summarized  in  Section  3.1.2.]  of  this  survey.  Wlien  these  same 
studies  were  performed  in  an  HLDCB,  the  trends  were  similar  to  those  of 
Section  3.1.2. 1 but  more  pronounced.  Effects  reported  Include; 

(a)  marked  rise  in  the  crack  pressure,  (b)  large  increase  in  propagation 
rate  as  crack  width  decreases,  and  (c)  complete  disappearance  of  the 
rising  (small  crack  width)  )jortion  of  the  curves  stiown  in  Figure  3,  even 
at  low  pressure.  Conclusion  (1)  of  Section  3.  1.2.1  (ignition  rate 
Increases  with  channel  Length  at  constant  lenglh-to-width  ratio)  and 
effect  (b)  noted  above  could  be  reconciled  If  it  were  known  that  the 
HLDCB  tests  had  been  run  at  constant  crack  length. 
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Bobolev  (1965a)  studii'd  the  flashdowii  into  a bJiiid  poro  in 

the  HLDCB.  The  results  are  preseated  in  Section  1.1. 1.2  where  the  blind 
and  open  pore  samples  burned  in  LLDCB  and  HLDCB  re.sults  are  compared. 

3 . 1 . 2 . 3 Blind  Single  Defect  in  AtmosphcT ic  Tests 

In  the  introduction  it  was  emphasized  that  convective  combustion  of 
porous  charges  and  single  defects  .irisoL  from  a pre.ssure  differential 
between  the  bomb  and  defect  pressures.  This  differential  consists  of 
two  components:  (1)  that  imposed  by  the  experimental  conditions,  and 

(2)  that  generated  by  combustion.  In  previous  sections  of  the  survey, 
the  first  contribution  has  always  outweighed  the  second.  An  .additional 
artifact  of  the  high  pressure  experiments  conducted  in  the  LLDCB  has 
been  the  dilution  effect  ol  pressurizing  gases  in  the  pores  of  the  test 
sample.  Elimination  of  this  problem  by  u.sing  an  "embedded"  charge  has 
only  emphasized  the  experimentally  created  pressure  difference. 

Godai  (1970)  used  two  defect  geometries  in  the  atmospheric  pres.suru 
blind  defect  pore  experiments:  (1)  .a  flat  crack,  5 mm  wide  by  25  mm 
long,  with  '.,’Acing  between  propellant  surfaces  van,'ing  from  0.05  to  0.5 
mm,  and  (2^  round  hole,  0.5  to  2.0  mm  in  diameter,  drilled  longitudinally 
through  a propellant  sample  5 x 5 x 25  mm.  The  propellant  formulation 
contained  75  or  80%  trimodal  AP  and  25  or  20%  polyester  binder.  Effects 
studied  included  propellant  temperature  and  pro]>ellant  formulation. 

The  effect  of  propellant  temperature  was  investigated  for  the  flat 
crack  at  -25°C  and  25°C  with  results  shown  in  Table  2.  Although  the 
trend  is  as  expected,  it  is  questionable  whether  the  precision  of  tlie 
reported  results  is  consistent  with  (a)  the  data  scatter,  admitted  by 
Godai  (1970),  and  (b)  tlic  effect  of  temperature  on  such  small  dimensions 
The  second  question  would  be  resolved  if  tlie  temperature  were  known  at 
which  the  gap  measurements  wore  made. 


TABLE  2.  Tliresliold  Conditions  for  Propag.it  i on  of 
Combustion  into  a Blind  Flat  Crack  at  Atmospheric  Pressure. 


Temperature,  °C 

Conductive  rate. 

-25 

0.77 

25 

0.  88 

mm/ s 

Threshold  ci 

0. 

0. 

widtli 


r.mi 


The  study  of  propel  l-ant  i orinulfitlon  v.irlables  included:  (1)  oxi- 

dizer particle  size  and  pei'centage,  (2)  addition  of  i-arhon  black,  and 
(3)  sub.sti  tut  ion  of  polybiiLadieiie  liimhr.  Tin-  c.'ffcct  of  oxidizer  pai-ti- 
cle  size  and  percentage  on  burning  rate  and  critical  cr.ick  gap  are  shown 
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in  Figure  6.  While  tlie  results  again  shoiv  the  expected  decrease  in 
threshold  crack  width  with  increasing  regression  rate,  it  is  clear  that 
burning  rate  alone  does  not  control  the  critical  conditions.  Part  of 
the  increase  in  threshold  crack  width  with  particle  size  may  be  attri- 
buted to  the  increased  roughness  of  the  channel  wall  at  larger  particle 
size.  This  would  have  two  effects:  (1)  increased  flow  resistance,  and 
(2)  either  .in  increase  or  decrease  in  channel  w.all  igni  tability,  depend- 
ing on  the  relation  between  particle  size  and  thermal  wave  thickness. 

Tlie  addition  of  carbon  black  produced  no  detect.ible  effect  on  the  rela- 
tionships shown  in  Figure  6.  It  was  concluded  that  neither  radiation 
from  carbon  particles  In  the  gas  phase  nor  increased  surface  absorptiv- 
ity contributes  to  flashdown.  Substitution  of  poiybutadieiie  binder  for 
polyester  resulted  in  lowering  of  the  threshold  crack  width,  which  is 
also  consistent  with  the  increased  burning  rate. 
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FIGURE  6.  Effect  of  Oxidizer  Particle 
Size  and  Percent.ago  on  Threshold  Crack 
Width  (Closed)  at  1 atm  (0.1  MPa).  [Data 
taken  from  Godai  (1970),  Figure  6.] 
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When  Godai  conducted  tests  on  the  blind  raund  liole,  flashdown  d i tl 
not  occur  tor  any  proi)e  1 1 an  t (ol  those  studied  in  tile  t Lai  crack),  even 
when  10  mm  holes  were  used.  The  ex[)lanaLion  ofti.Ted  was  the  possibility 
of  three  dimensional  recirculation  flows  in  the  1 lat  crack  which  tould 
not  take  place  in  the  round  tiole.  Considering  that  the  lar|',est  dimension 
in  the  flat  crack  was  only  5 imii,  one  would  expi'cl  that  recirculation 
effects  should  have  also  been  present  in  the  10  mm  round  liole. 

3.1.1  bxpe  r iiiieuts  Us  i ny  1 1^' i^je  n diiiy;Le  tore  JieoiiK^t  rv  (Macaron  L ) 

In  tills  configuration,  Finure  U-,  the  jiore  of  the  test  sample  is 
open  at  hotli  ends  thereby^  providing  a flow  channe  1 with  no  stagnation 
region . 

3 . 1 . 3 . 1 Macaroni  Sample  iii^  jy'"-'.  LLDCIi 

Althougii  initially  the  pore  ol  this  sample  is  filled  with  cold 
inert  gas,  there  Is  a vent,  other  than  the  entrance,  tliat  these  gases 
can  exit  through  when  the  eomhustion  ])eiU't  rates  the  pore'.  In  this 
configuration  the  sole  stimulus  causing  the  combustion  to  penetrate  into 
Che  pore  is  that  generated  hy  combustion. 

Margolin  and  Margulis  (1969)  carried  out  experiments  on  nitro- 
glycerin powder  in  several  e.onf Igurations  (Figure  7).  Fore  size  (1  mm 
diameter)  and  pressure  (50  kg/cm^)  were  chosen  by  Margolin  and  Margulis 

to  give  a value  of  An  (.see  Section  3.2.2.  1 for  discussion  of  Andreev 

number  criterion)  which  was  10-20  times  that  required  for  iiropagation  ot 
combustion  into  an  o|>en  ended  eiianneJ.  Results  for  the  condition:' 
depicted  in  Figure  7 were: 

a.  No  propagation 

b.  Propagation  only  as  tar  as  the  side  hole 

c.  initial  proi).ig.ti  1 on  into  both  ends,  finally  di:veioping  in 

only  one  side 

d.  Propagation  after  burnthrough  of  the  blind  end. 


FIGIIRK  7.  Pore  Con f i gu ra L ions  Used 
to  Study  Effect  of  Pore  Ceoiiietry  on 
Convective  liurning.  (Data  taken  from 
Margolin  and  M.argulis  (1969),  I’Lg- 
urc  1 . ] 
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Additional  cnperimLiilal  results  were  reported  by  Bobolev  (1965a) 
comparing  propagation  rates  into  blind  and  open  flat  cracks  ot  RDX 
(length  30  min,  particle  size  200  nm,  charge  density  1.2  g/cm^,  craik 
width  0.1  mm).  The  data  W'ln'  litted  to  tlie  toilowing  empirical  relation- 
ship.s : 

Open  crack 

...1.55 

V = 0. la  p (3) 

ih 

Closed  crack: 

= 0.075  p^ ■ (4) 

where  is  ignition  rate  in  cm/s  and  p is  in  atmospheres.  The  eiiect 

of  crack  length  and  width  was  not  studied;  from  tne  available  informa- 
tion it  is  seen  chat:  (1)  closing  the  crack  halves  tlie  late  ot  propa- 
gation, but  (2)  the  sensitivity  of  rate  to  pre.ssure  is  not  changed, 

Belyaev  (1969)  conducted  tests  on  secondary  explosives,  NO  powder, 
and  AP  mixtures  in  closed  and  open  flat  cracks  of  varying  widths. 

Belyaev  chose  conditions  to  insure  penetration  of  combu.scion;  however, 
the  value  ot  (see  discussion  in  section  3.2.2.  1 ri'Jating  to  Fispire 

31)  was  not  specified  so  time  the  margin  by  which  conditions  exceeded 
those  for  penetration  was  not  known.  Qualitative  comparison  of  ignition 
rates  for  open  and  closed  channels  is  shown  in  Figure  8.  Tlie  ignition 
velocity  and  maximum  crack  pressure  are  shown  in  Figure  9 as  functions 
of  crack  width.  No  propellant  formulation  was  given  for  cither  figure, 
nor  was  the  effect  of  Icngtli  on  ignition  rate  studied.  Experimental 
results  (Figure  10)  (still  on  unknown  propellant)  indicate  the  influence 
of  length  and  1 ength-to-width  rati<)  of  a closed  pore  on  ma.ximum  pore 
pressure' . A simp  Lifted  tlu'oretical  analysis  gives  trends  which  agree 
qualitatively  with  those  shown  in  Figure  10. 

Onset  and  propagation  of  eombustion  into  open  mided  pores  h.ivL-  Ih-l'ii 
mentioned  previously.  Qualitatively  it  was  noted  earlieT  [Margolin  and 
Margulis  (1969)]  that  spontaneous  pL-netration  into  an  eipi'n  pore  oecui's 
mucli  more  readily  than  for  a closed  pore,  the  threshold  Andn-ev  numtn'r 
(see  Section  3.2.2.])  being,  10-20  times  as  grL'.-it  for  the  closi.'d  i>ore. 

No  length  effect  was  rL’]iorted  by  Margolin  and  Margulis  (1969),  howi'viT, 
limited  data  [Bobolev  (1965.a)]  indicate  (see  Figure  11)  that  not  only  is 
the  threshold  pressure  lower  for  an  oj^en  channel  but  that  it  also  di'- 
creases  with  length.  DiffL^reilt  results  were  obtaiiiL'd  lor  rising  pii's- 
sures  as  also  shown  in  Figure  11. 

Ignition  jiropagation  rates  into  open  channe  I s of  KDX  wen-  studied 
and  rL'portL'd  | Bobolev  ( 1965a)]  (.sih'  Kq.  (3))  to  be  double  tliL-  value  tor 
closed  channels.  Similar  (jua  1 i t.i  t i ve  results  weri'  obtained  by  lU'lyacv 
( 1969). 


20 


[>ISlANk(  IC’.NlIfl)  f.NL) 


flGURt  8.  Quai i tati vt?  Comparison  of 
Convective  Burning  in  Ope.n  and  Closed 
Poren.  [Data  taken  from  Bclyaov  (196'1), 
Figure  3.  ] 


FlGl'RE  9.  Effect  of  Pore  Size  on  Maximum 
Channel  Pressure  and  Ignition  Patc‘.  (Un- 
known propellant  and  cliannel  lengLti.) 
[Data  taken  from  ilelyaov  (1969),  tabu- 
lated data . ] 
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Payii^  (1969)  ran  tests  cn  an  aluminized  AP  composite  at  700  psi^ 
(4.93  MPa)  in  nitrogen  using  cracks  constructed  of  slabs,  0.25  in.  long 
by  0.125  in.  wide  (6.4  x 3.2  inm) , and  having  spacings  rangii;g  from 
0.002-0.028  in.  (0.05  x 0.7  ram).  The  flame  entered  the  crack  in  all 
tests.  Time  required  for  combustion  to  enter  the  crack.s  Increased  with 
crack  width  while  ignition  rate  increased  as  crack  width  decreased. 

This  latter  relationship,  shown  in  Figure  12,  is  qualitatively  similar 
to  that  found  by  Belyaev  (1973)  for  clo..ed  cracks  (see  Figures  3 and  4). 


CPACK  win  in.  m.ri 


I I A I CHAl.  K 


FIGURE  12.  Effect  of  Open  Crack  Width  on 
Convective  Velocity  for  Composite  Propel- 
lant at  700  pslg  (4.93  MPa)  in  Nitrogen. 

(Data  taken  from  Payne  (1969),  Figure  4.] 

Prentice  (1962  and  1977)  has  studied  both  blind  and  open  single 
pores  in  the  LLDCB.  In  1962  he  studied  nitrocellulose-petrln  propel- 
lants. He  found  that  he  could  not  get  the  samples  having  a blind,  1/16- 
inch  (1.6  mm)  pore  to  fla.sh  at  any  pressure  up  to  550  psig  (3.9  MPa) 

(the  limit  of  his  apparatus).  Idien  he  burned  a modified  (pore  ignition 
by  burnthrough  of  a propellant  web  at  ignited  end  of  sample  instead  of 
by  hot  wire)  open-open  sample  (see  Figure  13)  the  data  pre.sented  in 
Figure  13  were  obtained.  The.se  data  show  flashdown  at  approximately  0.3 
MPa  bomb  pressure  for  the  open  system.  This  dramatically  show.s  that  it 
is  much  easier  to  flashdown  into  the  open-open  pore  than  it  i.s  to  fla.sh- 
dowu  into  the  blind  pore.  Prentice  also  found  that  f ashdoisni  will  occur 
at  pore  diameters  less  than  tlu;  IJame  stand-off  distance.  Additives  to 
the  nitrocellulose  propellant  significantly  altered  the  flaslidowu: 
catalyzed  propellants  f la.sb.ed  more  readily,  while  cooler  propellant.s 
(1,3,5%  paraformaldehyde  addition)  were  less  likely  to  flaslidown. 
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FIGURE  ]3.  Flashdown  Conditions  for  a 
Round  Open  Ended  Hole  in  NC/Petrin  Pro- 
pellant. [Data  taken  from  PieiiLioe 
(1962)  , Figure  4. ] 

Prentice  extended  his  earlier  work  to  opaque  composito  propellants 
(Prentice,  1977).  Both  blind  and  open  pores  using  previously  mentioned 
modified  ignition  system  were  studied  using  the  LLDCB.  Flashdown  was 
determined  through  the  use  of  Vibration  Response  .Spectroscopy  (VRS)  . 
Prentice  found  Chat  flashdown  was  facilitated  by  the  open-ended  pore  as 
opposed  to  the  blind  pore.  Figure  14.  Several  propellants  were  tested 
and  the  data  arc  presented  in  Figure  15.  Prentice  also  investigated 
double-base  propellants  having  a mesa  burning  rate  curve.  The  burning 
rate  curve  and  flashdown  curve  are  shown  superimposed  in  Figure  16. 

From  these  two  curves  Prentice  lias  concluded  that  those  factors  which 
affect  the  bunting  rate  also  affect  tlie  onset  of  flashdown. 

This  configuration  was  al.so  considered  In  a comparative  study  by 
Bobolev  (1.965a)  and  is  discussed  in  Section  3.  1.3.2. 

3 . 1 . 3 . 2 Experiments  Using  the  Open-Oiien  Single  Pore  r.eometry  (Macaroni) 
in  the  HL'dC^  ' ~ 

With  the  exception  of  gun  propellant  studies  which  used  perforaLed 
grains  burned  in  an  HLDCB,  very  little  work  has  been  iierformed  in  this 
area.  As  mentioned  previously,  Bobolev  (1965a)  compared  resul  t.s  for 
sampie.s  having  blind  and  open  pores  burned  in  both  LLDCB  and  HLDCB. 

This  work  presents  brief  experimental  results  of  onset  of  burning  in 
narrow  slits  made  from  RDX  and  Incite  (0.1  min  spacing).  Both  ri.sitig  and 
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FIGURE  16.  Burn  Hate  and  Flashdown  Condir.lon.s 
for  M-5  Propellant.  (Data  taken  from  Prentice 
(1977),  Figures  31  and  33.] 


constant  pressure  for  blind  and  open  slots  of  two  lengths  are  reported, 
with  the  results  shown  graphLcallj'  in  Figure  11.  General  conclusions 
are  that,  for  rising  pressure,  combustion  penetrates  a closed  channel 
more  readily  than  an  open  one;  on  the  other  hand,  for  constant  pressure, 
penetration  is  facilitated  by  an  open  channel.  In  addition,  lengthening 
the  channel  makes  penetration  easier  except  for  an  open  clianncl  witli 
constant  pressure  environment.  Since  the  data  given  are  so  sparse  and 
appear  to  be  somewhat  contradictory,  with  no  principles  elucidated  to 
explain  the  results,  the  value  of  this  paper  is  that  it  raises  questions 
which  should  be  answered  by  future  experiments  and  analysi.s.  It  would 
be  especially  valuable  to  repeat  the  tests  for  otiicr  propellants  and 
rates  of  pressure  rise  (.none  were  specified  by  Bobolev  (1965a)). 

3 . 1 . 3 . 3 Experiments  Using  the  Open-Open  Single  Pore  Geometry  (Hac:aronl) 
Burned  at  1 Atmosphere 


Whereas  the  blind  pore  experiments  performed  at  1 atm  (0. 1 HPn) 
were  run  mainly  on  flat  cracks,  nearly  all  the  open  pore  tests  utilized 
the  round  hole  (Godai  (1970)].  Effect.s  studied  were  (1)  oxidizer  parti- 
cle size  and  percentage  (AP/pol yes  ter  formulation),  and  (2)  carbon  and 
copper  chromite  F additives  (AP/poiybutadiene  formulation).  The  in- 
fluence of  AP  particle  size  and  percentage  is  shown  in  Figure  17.  The 
particle  size  dependence,  noted  for  tlie  blind  crack,  is  still  apparent; 
however,  tlie  percentage  of  AP  has  practically  no  effect.  The  effect  of 
addition  of  aluminum  on  burning  rate  and  critical  diameter  is  .shown  in 
Figure  18.  The  interesting  result  here  is  that  conductive  rate  is 
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FIGURE  17.  Effect  of  Oxidizer  Particle 
Size  and  Percentage  on  Threshold  Pore 
Diameter  (Open)  at  1 atm  (0.1  MPa). 

[Data  taken  from  Godai  (1970),  I’igurc  7.) 


ALUMINUM.  V 

FIGURE  18.  Effect  of  Aluminum  on  Conductive  Rati- 
and  Thre.shold  Pore  Diameter  for  Open  Round  Pore. 
[Data  taken  from  Godai  ( 1970),  Figuri'  8.] 
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virtually  uiiaf  ferted  while  flaslidown  diameter  dopeitds  cm  both  oxidizer 
and  aluminum  content.  Again,  there  i.s  clear  evidence  that  conductive 
rate  alone  does  not  determine  flashdown  condition.  The  ri'sulL  of  .adding 
the  burning  rate  catalyst,  copper  chromite  F,  is  .shown  in  Figure  19  and 
suggests  that  conductive  rate  is  at  least  one  ot  the  contributing  i ac- 
tors in  flashdown. 


0 i I 1 1 I 

0 O.i  1 1.^  ? 

COPPER  (.HPOMITE.  % 


ROUNO  HtILfc. 


FIGURE  19.  El  feet  of  Copper  Chromite  on 
Conductive  Rate  and  Threshold  Port>  Ui.ameter 
for  Open  Round  Pore.  (Data  taken  from 
Godai  (1970),  Figure  9.] 

No  details  were  given  of  the  results  of  flashdown  into  open  flat 
cracks.  It  was  noted  that  (1)  the  flame  velocity  accelerated  more  tlian 
for  closed  flat  cracks,  and  (2)  the  threshold  crack  width  was  approxi- 
mately 107,  less  for  the  open-ended  case.  This  latter  ob.servat ion  pro- 
vides additional  evidence  for  a reeirculatiug  flow  in  Llie  elo.sed,  flat 
crack. 


3.  1.3.4  Experiments  Using  the  Opc-n-Open  Single  Pore  Geometry  in  a 
Special  Vented  LLUCB 

The  development  of  combustion  in  a single  pore  or  in  a porous  bed 
first  requires  the  penetration  of  liot  combustion  gases  evolved  in  the 
normal  burning.  Secondly,  these  1 lowing  gases  must  heat  the  pore  wall 
to  a condition,  usually  called  ignition,  wliieh  leads  to  self-sustained 
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combustion.  Because  tliere  is  a tinu'  delay  between  penetration  and 
ignition,  the  entering  gases  flow  ahead  of  the  ignition  iront. 

Krasnov  (1970)  devised  an  experiment  enabling  the  measurement  of 
both  the  velocity  of  inflow  and  velocity  of  propagation  of  the  ignition 
front.  A sample  of  balllstite  powder  (double-base),  containing  a single 
cylindrical  pore,  40  mm  long  and  of  varying  diameters,  was  ignited  in  a 
constant  pressure  burner.  The  unignited  end  of  the  pore  was  connected 
to  the  atmosphere  through  a fi.xture  containing  a solenoid  valve.  Tvyo 
different  approaches  were  used.  In  one,  the  pore  was  initially  filled 
with  a liquid  such  as  water  or  alcohol;  upon  ignition,  the  solenoid 
valve  was  opened  and  the  liquid  was  expelled  at  a rate  depending  on  the 
valve  area.  Both  the  liquid  meniscus  and  ignition  front  were  recorded 
by  cinephotography  through  the  relatively  tr-ansparent  proi^eliant.  In 
the  other  approach,  a flow  of  nitrogen  was  established  through  the  pore 
prior  to  ignition  by  adjusting  the  valve.  The  two  methods  were  stated 
to  give  identical  results,  thereby  presumably  dispelling  any  doubts 
regarding  inertial  or  cooling  effects  of  the  liquid. 

Trends  of  experimental  results  are  shown  in  Figure  20  for  an  un- 
specified pressure.  The  ordinate  scale  shows  the  ratio  of  gas  flow 
velocity  to  ignition  front  velocity  for  pore  sizes  ranging  from  I to  4 
mm.  Both  velocities  were  tak(!n  from  tlie  cinephotographs  vyith  v^^.  being 
the  velocity  of  the  liquid  meniscus.  In  general,  the  ignition  lag 
increases  with  increasing  flow  velocity  and  decreasing  poi'e  size.  For 
the  conditions  of  the  ex])eriment,  the  ignition  lag  is  independent  of  pore 
size  above  2 mm  and  also  of  flow  velocity  above  a critical  value  which 
increases  with  pore  diameter.  The  leveling  off  of  ignition  lag  with 
increasing  pore  size  is  shoiyn  in  Figure  21,  which  is  cross-plotted  from 
Figure  20  for  a flow  velocity  of  30  cm/s.  It  is  possible  to  account  for 
the  effect  of  pore  size  if  it  is  noted  tliat  the  velocity  of  the  ignition 
front  increases  for  pores  larger  than  2 mm  diameter.  Auxiliary  experi- 
mental results  showing  this  effect  are  shown  in  Figure  22  for  a pressure 
of  40  atm  (4  MPa).  Even  though  Figure  22  represents  experiments  relevant 
to  a blind  pore,  the  data  are  probably  applicable  here  because,  as  indi- 
cated by  the  consistent  rise  in  ignition  velocity  with  pore  diameter,  the 
dead  end  effect  is  negligible  in  this  (.-xper  iment . 

The  underlying  reason  for  the  effect  of  diameter,  not  mentioned  by 
Krasnov  (1970),  is  connet:ted  with  the  relation  between  normal  f 1 aiiR' 
standoff  distance  and  pore  size.  At  constant  pressure,  this  distance  is 
also  constant,  so  that  larger  pores  permit  the  entry  of  the  higher  Lemper 
ature  gases  associati-d  wtli  complete  comljustion.  It  would  he  in  Li' res  t i ng 
to  repeat  the  experiment  for  a composite  propellant,  for  whicli  the  gas 
phase  reactions  are  completed  closer  to  tlie  surface  than  for  balllstite. 
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FIGURE  20.  Ratio  of  Velocities  of  Gas 
Flow  (Vg(.)  and  Ignition  Front  (vig)  as 
a Function  of  Gas  Flow  Velocity  and 
Pore  Diameter,  dp  (Shown  on  Curves)  . 
[Data  taken  from  Krasnov  (1970) , 

Figure  4. ] 


FIGURE  21.  Ratio  of  Velocities  of  Gas 
Flow  (Vgj.)  and  Ignition  Front  (tig)  as 
a Function  of  Channel  Diameter,  Vg(,  =• 
30  cm/s.  [Data  taken  from  Krasnov 
(1970) , Figure  5. ] 


NWC  TP  6007 


FIGURE  22.  Effect  of  Channel  Diameter  on 
Rate  of  Flow  Into  Blind  Channel  at  40  atm 
(4  MPa).  [Data  taken  from  Krasnov  (1970), 
Figure  7 . ] 


This  completes  the  survey  of  literature  dealing  with  early  stages 
of  convective  burning  in  single  defects  in  LLDCB,  HLDCB,  and  at  atmos- 
pheric pressure.  Pressure  effects,  but  no  quantitative  effects  of  the 
rate  of  pressure  rise,  have  been  noted.  Thi-  next  section  deals  with 
results  of  experiments  conducted  on  porous  beds. 

3.2  EXPERIMENTS  USING  THE  POROUS  BED  GllOHETRY 

While  samples  containing  single  pores  may  be  char.ac ter izi-d  by 
rather  precisely  measured  sizes  such  as  lengtli  and  cross-set:tional  area, 
the  porous  samples  must  be  characterized  by  averaged  quantities  sucli  as 
porosity,  average  particle  size,  specific  surface  area  and  permeability 
to  flow.  Porous  beds  and  single  pore  samples  differ  by  degree  ratlier 
than  kind.  But  since  they  do  represent  opposite  ends  of  the  increased 
surface  area  spectrum,  tliey  will  be  discussed  .separately.  Tlie  samples 
having  single  pores  or  channels  have  already  been  discussed.  The  en- 
suing discussion  will  deal  with  porous  beds. 

3.2.1  Embedded  Porous  Bed  Burned  in  the  LLDCB 

Both  onset  and  development  of  convective  burning  have  been  .sLudii’d 
for  porous  beds  initially  sealed  against  the  bomb  pressure.  No  detailed 
presentation  was  made  of  tne  criteria  lor  the  onset  phase;  instead,  it 
was  simply  stated  [Bely.acv  (1966)]  that  results  were  similar  to  those 
obtained  using  the  HLDCB  witli  a closed  porous  bed  (see  Sectiou  3. 2. 2. 2). 
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The  implicntion  is  that  tho  onset  i>hase  ot  convective  burnini',  di'iieiuls 
only  on  a preasnre  differenci-  between  the  bomb  and  porous  bed  iiressurcs, 
whether  the  difference  arises  from  an  increasing  bomb  pressure  or  t rom 
an  initially  high  bomb  pressure.  It  would  be  expected  that  the  rate  of 
pressure  rise  for  the  HLDCli  would  iiifluenee  the  results,  but  the  eftect 
appears  not  to  have  been  studied. 

The  rate  of  projiagation  of  comhuslioii  into  an  embedded  porous  bed 
burned  in  an  LLDCB  was  the  subject  of  two  papers  iDubovitskii  (1974a) 
and  Frolov  ( 1972)].  In  addition,  a third  reference  [Uubovitskii  (1974b)  ] 
apparently  a different  traiislation  of  Diibovi  tski  i (1974a)  was  included 
in  the  survey  in  order  to  resolve  certain  translation  discrepancies. 
Initial  conditions  of  the  expiT imeiit  were  chosen  to  t.'nsure  the  onset  ol 
convective  burning.  Effects  investigated  were:  propellant  type,  flame 
temperature,  porosity,  (larticle  size,  .and  sample  length.  The  position 
and  temperature  of  tlu’  advancing  combustion  front  were  measured  oiitiral- 
ly  and  by  means  of  thin  wire  tlierraocou|) les  embeddc’d  at  right  angles  to 
the  direction  of  propagation. 

Experiments  were  conducted  at  an  initial  bomb  pressure  of  70  atm  (7 
M?a)  on  an  RDX/.<\P  composition  (Dubovitskii  (1974a)]  ha^’inp,  ?07  porosity. 
The  velocity  of  the  advancing  combustion  iront  is  shown  as  a lunction  ot 
its  position  in  Figure  23.  Time  also  varic's  along  eacli  curve  hut  is  not 
indicated  in  the  original  reference.  IVo  A?  particle  size  ranges  (100- 
125  urn  and  200-250  urn)  were  tested  but  no  mention  was  made  of  the  size 
range  used  in  obtaining  the  data  for  Figure  23;  it  was  stated  that 
propagation  velocity  increases  with  increasing  particle  size.  Flame 
temperature  was  varied  by  changing  proportion  of  AP  and  RDX.  Test 
results  can  be  summarized  In  the  following  statements: 

1.  Combustion  propagation  velocity  is  greater  for  longer  charges, 
the  effect  becoming  more  pronounced  at  greater  distances  from  the  ig- 
nited end  of  the  charge. 

2.  Combustion  propagation  velocity  decreases  as  the  end  of  the 
sample  is  approached. 

3.  Combustion  propagation  velocity  in  Region  1 is  greater  for 
higher  flame  temperature,  increasing  with  distance  for  low  flame  temper- 
ature, and  decreasing  with  increasing  distance  for  high  flame  temper- 
ature . 


The  first  two  effects  are  a re.sult  of  tin-  counter  pressure  produced 
as  the  combustion  front  approactie.s  tlie  closed  end.  The  increase  in 
propagation  velocity  with  combustion  temperature  is  readily  explained  by 
the  greater  ease  of  Ignition  of  the  jiore  walls  by  the  hotter  gases.  The 
reason  for  the  increase  in  velocity  with  distance  for  low  gas  tempera- 
ture compared  with  a decrease  lor  ii  i gli  );aH  temperaturi'  is  not  (xijlained 
in  a satisfactory  manner.  It  ap])ears,  however,  to  be  the  result  oi  an 
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FIGURE  23.  Velocity  of  Convective 
Combustion  Front  vs.  Position  for 
RDX/AP  Porous  CiiarRc.s.  [Data  taken 
fron  Dubovitskii  (1974a),  Figure  1.] 


interaction  between  gas  flow  rate  and  ignition  time  as  gas  temperature 
increases.  In  examining  Figure  23,  one  notes  three  distinct  regions  of 
convective  propagation.  The  first  region,  adjacent  to  the  ignitea  end, 
is  a region  of  decreasing  velocity  in  which  gas  flow  and  heat  transfer 
are  the  principal  effects.  Chemistry  and  ignition  have  not  yet  bccotiK- 
important;  extent  of  Region  1 appears  to  be  insensitive  to  sample  length 
(for  lengths  cm).  Tne  second  region  is  one  of  increasing  velocity 
following  ignition.  The  extent  of  Region  11  and  acceleration  of  gases 
both  increased  with  sampli’  lengtit,  v.;Ij1c1i  determines  how  soon  the  pres- 
ence of  the  closed  end  affects  the  progress  ol  convective  burning. 

Finally,  in  Region  111,  the  gases  slow  down  owing  to  the  approach  to  the 
closed  end  and  the  development  of  counter  pressure. 

In  addition  to  studies  of  the  velocity  of  propagation  of  the  com- 
bustion front,  measurements  were  made  of  the  maximum  pore  ana  l>omb 
pressures  as  a function  of  sample  length  and  AP  particle  size  i Duhov i tski i 
(1974a)].  Results  are  ])lotted  in  Figuia'  24,  indicating  the  rather 
significant  increase  in  maximum  pori-  pri.‘ssure  and  moderate  risi'  in 
maximum  bomb  pressure  with  s-amplo  length.  The  pore  jiressure  decreases 
with  increasing  particle  size  bi'causc,  for  the  same  porosity,  tlie  ri'sisL- 
ance  to  flow  Is  less  for  larger  particle  size. 
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FIGURE  24.  Maxiniuiii  I’ort.'  and 
H(imb  ProsBurt's  I t'r  Cnnvoot  ivo 
Burning  in  RUX/AP  Porous  Cliarges. 
(Data  taken  from  OuboviLskii 
(1974a).  Tabio  1.] 


Similar  tests  wore  conducted  and  reported  [Frolov  (1972)]  for  a 
stoichiometric  mixture  of  AP  and  polystyrene  (particle  size  ranges  bO- 
100  nm  or  100-130  nm)  having  a flame  temperature  of  2800°K.  The  chargi-s 
were  all  10  mm  in  diameter,  20  and  307’  porosity,  and  were  burned  nr  an 
initial  bomb  pressure  of  50-70  atm  (5-7  MPa)  with  pressure  rises  of  5- 
lOZ.  Velocity  of  the  gas  front  as  a function  of  position  is  plotted  in 
Figure  25.  The  general  shape  of  the  curves  is  similar  to  those  sliown  in 
Figure  23,  in  addition  to  which  it  is  seen  that  pio[agaLion  velocity  is 
greater  at  high  porosities  but  that  the  flow  decelerates  more  rapidly  in 
Region  I for  low  porosities.  Comparison  of  Figures  23  and  25  fiugge.st.s 
that  the  propellant  chemistry,  especially  ignitability  which  de])cnds  on 
kinetics,  has  a significant  influence  on  the  e.xtent  of  Region  1. 

3.2.2  Blind  Porous  Beds 

3. 2. 2.1  Blind  Porous  Bed  .Burned  in  the  1.1. DCB.  This  section  sum- 
marizes results  obtained  by  burning  porous  beds  in  a Crawford  type  bomb 
(LbDCB)  in  which  the  high  r.ates  of  pressure  rise  associated  with  the 
manometric  bomb  (HLDCB)  are  absent.  Since  tin-  defect  is  exposed  direct- 
ly to  the  igniter,  special  efforts  were  usually  made  to  provide  a "soft" 
Ignition,  thereby  avoiding  the  effect  of  a rapid  pressure  line  during 
the  ignition  transient.  The  degree  of  succes.s  ot  these  attempL.s  cannot 
be  assessed  because  the  cili-ct.s  of  igniter  size  and  brisanee  were  not 
determined. 
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FIGURE  25.  VoLocity  ol'  G^>nv^'^T  i vl'  Combust  iun 
Front  vs.  Position  S toifli  i onu’ I r i r AP/Po  Ly.styronL! 
Porous  Cliarjios.  [Data  l^)kon  t rom  Frolov  0972), 
Fij;uro  2.  ] 


AJ  thouslt  the  oxperimont.il  results  to  In-  summirlxod  in  this  section 
relate  to  geometricall v blind  porous  beds,  no  studies  wito  reported  ol 
the  effect  of  the  length  of  the  test  sample.  Ui'in-c'  it  is  not  possible 
Co  resolve  whether  or  not  the  samples  were  hydrodynami e.illy  blind.  In  a 
parallel  effort  involving  rising  pressure  (Section  3. 2. 2. 2)  the  length 
was  adjusted  so  that  its  effect  w,-is  minimized.  Perhaps  the  same  jiro- 
cedure  was  followed  for  the  Cr.awford  bomb  tests.  It  is  again  euiph.isi zed 
that  in  these  tests,  the  porous  bed  is  open  to  tlie  bomb  pressure  prior 
to  ignition  so  chat  the  pores  are  iilJed  with  higli  pressure,  et)ld,  inert 
gas. 


Of  tlie  references  found  ri'levant  to  onset  of  eonvt'ctive  burning  in 
blind  porou.s  beds,  one  presents  no  ex]ier imental  data;  the  statement  is 
made  that  threshold  pressures  were  10—15  times  as  gre.at  as  Lliose  ob- 
tained in  experiments  using  the  manometric  bomb  (Section  3. 2.2. 2)  oi‘  the 
Crawford  bomb  with  "embedded''  charges  (Si'ction  3.2.1).  It  similar  LUid 
effects  are.  assumed  in  all  three  series  ol  ti’sts  (all  were  geometrically 
blind)  then  the  greater  stability  to  onset  of  convective  burning  in  llie 
present  case  is  explainable  bv  tlu'  initial  presetice  of  the  cold,  inert, 
'jure  gases.  Tliese  cold  gasi-s,  at  bomb  pressure,  both  impeded  the  ini  low 
of  combustion  products  and  cooled  tin-  combustion  products  thar  did  i low 
in  so  that  pore  ignition  w.is  more  dil  ficulL. 

The  most  significant  result  ot  studies  of  onsi’t  of  combustion  in 
porous  beds  by  spontaneous  penetration  was  the  confirmation  of  ,i  concept 
whicli  has  been  referred  to  as  the  "Andreev"  criterion.  A diversion  to 
present  .!  brief  background  ol  tin-  Andreev  criterion  is  appropriate  lunc. 
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It  is  jihysic.illy  reasonable  to  say  tliat  combust  i on  of  a porous 
material  will  proceed  in  a normal  manner  as  long  as  the  surface  irregu- 
larities, resulting  from  the  porosity,  are  smaller  than  some  dimension 
reflecting  the  combustion  process.  Helyaev  proposed  this  idea  in  an 
unavailable  doctoral  thesis  in  1946.  He  suggested  using  the  ratio  of 
pore  size  to  standoff  distance  of  the  gaseous  reaction  zone  as  a criter- 
ion, probably  with  the  idea  of  application  to  a single  pore.  Tlie  idea 
was  quantified  by  Andreev  in  1946  in  another  un.available  paper  still 
applied  to  a single  jiore  and  using  the  solid  phase  Lhermai  layer  thick- 
ness as  the  combustion  dimension.  Tlie  first  evidence  linking  the  cri- 
terion to  a statistically  porous  bed  X\;as  presented  by  Margolin  (1961)  in 
a short  theoretical  note  witliout  experimental  rt'sulis.  The  current 
status  of  the  Andreev  number  criterion  is  as  follows: 

1.  Although  Initially  developed  for  single  pores,  the  Andreev 
number  criterion  has  not  been  extensively  tested  experi- 
mentally for  single  pores.  (Exceptions:  see  Prentice  (1962 
and  1977 J (Section  3. 1.3.1  of  this  survey)  for  single  pore  re- 
sults . ) 

2.  i'he  Andrtu'v  number  criterion  has  been  aj^plied  only  to  spontan- 
eous penetratioii  of  combustion  into  single  pores  and  porous 
beds,  under  experimental  conditions  similar  to  those  found  in 
the  LLUCB  with  pores  prepressurized  to  the  initial  bomb  pres- 
sure . 

3.  The  Andreev  number  criterion  has  been  applied  to  porous  beds 
by  use  of  an  equivalent  pore  diameter  based  on  permeability 
measurements . 

4.  There  have  been  no  compreljeiisive  reports  of  experiments  de- 
signed to  determine  tlie  quaiitiLative  effects  on  Andreev  number 
criterion  of  end  condition  (open  or  blind)  or  of  length  of 
single  pores  or  porous  beds. 

Quantitatively,  the  Andreev  criterion  takes  the  form 

pre  d/A  = An  > An.  (5) 

s s * 

ivhere 

p = bulk  density  of  solid  (mass/ total  volume) 

Cj^  = specilic  lieaL  of  .solid 

d = cliaracter ist ic  defect  size  (actual  or  hydraulic  diameter) 

Xg  = thermal  conductivity  of  solid 
An  = Andreev  number 
An^  =»  Llireshold  Andreev  number 
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Rctur'vin,^  lo  specific  papers  iii  Lhe  literature  survey,  one  finds 
two  references  dealing  with  spontaneous  penetration  of  combustion  into 
porous  beds.  The  first  [Margolin  and  Chulko  (J.966)]  represents  a com- 
pendium of  experimental  results  obtained  for  single  component  materials 
having  a narrow  particle  size  range  and  taken  from  many  sources. 
Margolin  and  Chuiko  (1966)  presented  Eq . (5)  in  the  alternate  form 


p^rd  (6)  I 

in  order  to  eliminate  the  specific  heat  and  thermal  conductivity,  which  S 

are  not  well  known,  but  which  are  stated  to  vary  over  narrow  limits.  It  | 

is  Important  to  note  that,  while  An  is  dimensionless  for  consistent  | 

units,  <(1  has  dimension.s  of  g/cm-s  in  the  units  chosen  (p;g/cm3,  r:cm/s,  j 

d:cm) . In  applying  Eq . (5)  or  (6)  to  a single  pore,  an  appropriate  fl 

value  of  d would  be  the  pore  diameter.  In  the  case  of  a porous  bed,  the  I 

mean  hydraulic  pore  d'  imeter  is  used.  In  an  article  by  Margolin  and  | 

Chulko  (1966),  where  particle  size  varic.s  over  a narrow  range  for  any  i 

given  test,  the  mean  hydraulic  diameter  is  computed  from  J 

d “ (2/3)D  (1-.S)/.S  (see  Belyaev  (1973),  (7)  j 

Section  4)  3 

j 

wherf  -5 

Dp  = particle  di.ameter  | 

6 = relative  density  (one  minus  porosity)  s 

Materials  included  in  the  study  were  TNT,  picric  acid,  nitrocellulose,  J 

PETN,  nitroglycerin  powder,  tetryi,  RDX,  ItMX,  AP,  and  mercury  fulminate.  1 

Particles  range  in  size  from  20  to  850  urn  with  porosities  of  from  25.5  I 

to  87.5%.  Threshold  values  of  mass  burning  rate,  (pr)jij  and  hydraulic  j 

diameter  d^  are  shown  in  Table  3 for  two  pressures,  1 and  10  MPa.  Also  I 

shown  for  each  pressure  is  the  ratio  of  critical  diameter  to  critical  4 

diameter  for  TNT.  This  number  serves  to  rank  materials  according  to  | 

their  propensity  to  develop  convective  burning,  or,  in  other  words,  to  ^ 

their  stability  of  normal  burning.  The  higher  the  ratio  just  defined,  ] 

I he  lower  the  stability.  It  should  be  noted  that  the  order  of  stability  1 

'hunges  with  pres.sure  owing  to  the  nonlinc.ir  dependence  of  conductive 
burning  rate  on  pressure;  altliough  this  method  of  ranking  is  u.seful,  it 
is  important  to  specify  tlie  pressure.  It  is  intei'esting  to  note  that 
both  .'J  and  IIMX,  co ingredients  of  many  modern  high  energy  [)roiie]  lants, 

exliibit  among  the  highest  tendencies  to  develop  convective  burning.  In  g 

addition,  It  has  been  noted  [Kraeutle  (March  1974)]  that  tlie  exothermic  | 

decomposition  temperature  of  HMX  is  significantly  lowered  in  compo.sl- 

tions  containing  AP.  1 
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TABLE  3.  Threshold  Conditions  for  Breakdown  of 
Conductive  Burning." 


Explosive 

= (Prd)^ 
mg/cm-s 

Threshold  pore 
diameter,  d^  pm 

Ratio  c 

^*'^^*TNT 

1 fU'a 

10  MPa 

1 MPa 

10  MPa 

TNT 

8.0 

850 

108 

1.00 

1.00 

Picric  acid 

8.  5 

625 

90 

1.36 

1.20 

Nitrocellulose 

9.8 

330 

100 

2.58 

1.08 

PETN 

11.1 

540 

55 

1.57 

1.96 

Nitroglycerin 

powder 

4.35 

120 

40 

7.08 

2.70 

Tetryl 

7. 1 

235 

48 

3.62 

2.25 

RDX 

6.9 

163 

21 

5.21 

5.  14 

HMX 

6.5 

140 

24 

6.07 

4.50 

AP 

3.7 

110 

20 

7.73 

5.  40 

Mercury 

fulminate 

11. A 

115 

4.3 

7.39 

9.39 

^Slightly  different  values  are  given  by  Belyaev  (J973)  (Table  6). 


Noting  the  relatively  narrow  range  of  values  of  Margolin  and 

Chuiito  (1966)  took  its  average  value  (9  x 10"^  g/em-s),  c'oinblned  it 
with  a typical  value  for  A^/cj.  for  the  condensed  phase  (1.5  x 10~3 
g/cra-s),  and  so  derived  the  threshold  Andreev  number  (An^  = 6).  Owing 
to  the  simplicity  of  the  concept  and  to  the  averaging  technique  used  in 
deriving  the  threshold  constant.  It  is  nut  surprising  that  deviations 
arc  found.  Indeed,  Margolin  and  Chuiko  (196b)  report  variations  in 
which  <J)j^  increases  by  a factor  of  two  for  very  smooth  surfaces  and 
decreases  by  a factor  of  one-half  for  very  fine  particles.  In  addition, 
the  effect  of  combustion  temperature  and  flame  zone  structure  liave  been 
found  important.  For  example,  TNT,  with  the  highest  value  ol  also 

has  the  lowest  combes t Ion  temperature  ol  tlie  reported  secondary  explo- 
sive. Certainly  th-  implicit  effects  of  particle  size  and  fl.ime  temper- 
ature on  the  quantities  wlilch  comprise  and  An^  are  not  comiilc  tely 
expressed  in  the  simplified  criteria.  The  authors  recognized  this 
limitation  when  they  recommended  determination  of  An^  (or  U ' new 

classe.s  of  materials  rather  tiian  indiscriminate  application  ot  An^^>6.  A 
controversial  and  unresolved  issue  is  raised  by  tlie  claim  that,  for  a 
pore  size  di.stribution , the  appi’opriati'  v.ilue  of  diameter  to  be  used  in 
computing  Andreev  number  should  be  based  on  maximum  pore  size.  Otlier 
presentations,  such  as  the  one  which  follows  [Bobolev  (1966)],  conclude 
lliat  the  average  port.'  size  i.s  more  relevant. 
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The  otlicr  reference  related  to  Andreev  number  [Bobolev  (1966)] 
presents  results  of  studies  on  RDX,  AP,  and  mixtures  consisting  of 
RDX/AP,  AP/ uro tropine , AP/sucrose,  potassium  perchlorate  (KP) / uro- 
troi)ine,  and  KP/sucrose.  Wliilc  Margolin  and  Chuiko  (L966)  dealt  with 
single  component  systems  of  uniform  particle  size,  Bobolev  (1966)  was 
concerned  with  establishing  convective  burning  criteria  for  mixtures  and 
multimodal  particle  size  distributions.  In  addition,  the  effect  of 
component  melting  is  included  in  the  final  four  mixtures  since  AP  and 
urotropine  do  not  mcit,*  while  KP  .and  sucrose  do. 

The  effect  of  particle  size  and  pressure  on  mass  burning  rate  of 
RDX  (theoretical  maximum  density  = 1.80  g/cm^)  is  shown  in  Figure  26,  in 
which  the  lower  curve  shows  the  normal  burning  rate.  The  stand.ard 
Russian  method,  given  by  Belyaev  (1973)  for  calculating  mass  burning 
rate,  is  to  multiply  the  observed  conductive  rate  by  samplt'  bulk  density. 
Since  sample  bulk  density  varied  only  slightly  (0.996-1.038  g/cm^)  in 
this  series  of  tests,  Figure  26  sh.ows  mainly  the  effect  of  particle 
size.  The  high  slope  curves  indicate  convective  burning  rates  with 
particle  size  increasing  toward  tlu-  left.  For  larger  particles,  the 
interconnecting  pores  are  larger,  thereby  allowing  penetration  of  high 
temperature  combustion  products  and  development  of  convective  burning 
(intersection  point  witii  lower  line;  at  lower  pressures.  Conditions  for 
onset  of  convective  burning  of  both  pure  RDX  and  pure  AP  are  shown  in 
Table  4 for  different  particle  size  ranges.  The  hydraulic  pore  diameter 
used  in  the  calculation  of  v,.,  was  computed  as  the  mass  weighted  h.arinonic 
mean  of  the  particle  size  distribution  (see  Eq . (9)  below).  The  ration- 

ale for  this  approach  is  not  explained,  nor  is  the  final  value  of  pore 
diameter  which  lies  outside  tiu'  particle  .size  range.  Acceptance  of  the 
results  shows  that  for  eacli  substance  there  is  a near  constant  value  cjf 
and  that  normal  conductive  burning  of  RDX  is  more  stable  than  Lh.aL 
of  AP, 


Additional  sLudie.s  were  c.onductc'd  on  samples  of  RDX  conL.nining 
several  ranges  of  particle  sizes  (bimodal  and  trimodnl).  Three  alterna- 
tive methods  wore  used  for  c:omputing  the  mean  hydraulic  pore  diameter  to 
be  used  in  evaluating  fl-^; 


d =^a.D  . 

1 / I 1 pi 

d.  = 1/d,,  = Va,/D  . 

2 1 pi 

- xflkTTT^ 


(8) 

(9) 

(10) 


where  the  a, 
1 


are  weight  fractions  of  partiiles  with  diameti^rs, 


D . 

pi 


:*r 

Hightower  ,ma  Price  (1967  and  1968)  .and  Bogg.s  and  Kr.aeutle  (1969) 
have  shown  that  /\1’  has  a thin  liquid  layer,  several  mierons  thiik  on  tlie 
burning  surface. 
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Tho  origin  of  Eq . (8)  and  (9)  i.s  not  provided;  the  pore  diameters 
dq  and  d2  represent  arithmetic  and  harmonic  mean  mass  averaged  particle 
diameters.  Neither  d^  nor  leads  to  a constant  computed  value  of 
However,  the  hydraulic  diameter  83,  determined  by  measurements  of  the 
gas  permeability  k,  does  yield  a near  constant  value  of  independent 

of  the  particle  size  distribution.  Figure  27  compares  the  6^  values  as 
calculated  from  dq,  82,  and  83  for  bimodal  RDX  containing  the  particle 
size  ranges  50-63  um  and  160-250  urn.  Addition  of  a larger  size  range 
(250-315  um) , to  give  a trimodal  sample,  does  not  significantly  change 
the  threshold  value  if  83,  based  on  permeability  measurements,  is 
used.  Bobolev,  and  others  (1966)  conclude  that  the  constant  value  of 
based  on  permeability  measurements  means  that  the  ]>oreH  of  average 
diameter,  being  the  most  numerous,  are  the  defects  that  influence  onset 
of  convective  burning.  Tliis  is  in  contrast  to  Margolin  and  Chuiko 
(1966),  in  which  pores  of  maximum  size  were  postulated  as  governing. 

The  final  studies  described  by  Bobolev  (1966)  concern  the  influenc.e 
of  composition  on  Test  series  were  run  to  examine  the  following: 

1.  Additivity  of  for  mixtures  containing  the  two  monopropel- 
lants AP  and  RDX. 

2.  Effect  of  fu.sibility  of  the  fuel  and  oxidizer  eompuueiiLs  of 
the  mixture. 


100:0  7i>;2b  bO:bO  0.100 

^^APIICLt  mass  ratio.  SO-63  iriicron/ 1 bO  P'aO  inicroii 


computed  from 

(1)  d (Eq.  (8)) 

(2)  dj  (Eq.  (9)) 

(3)  d (Eq.  (10)) 


FIGURE  27.  Thre.sliold  Quantity  for 
RDX  vs.  Mass  Ratio  of  Fine  to  Coarse 
Particles.  (D.3ta  taken  from  Bobolev 
(1966),  Fiq;ure  5.] 
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The  AP/RDX  results  are  sunmiarized  in  Table  5,  which  sliows  the 
effect  of  mixture  ratio  on  iji^  (calcuLateJ  Ltum  il_j)  for  a particle  size 
range  of  160-250  nm.  The  interesting  result  is  that  does  not  seem  to 
be  an  additive  function  of  the  ingredients  (tl)is  would  yield  a linear 
change  in  with  percent  AP)  but  changes  abruptly  at  that  composition 
where  the  volumetric  fraction  of  RDX  exceeds  50%.  No  explanation  of  the 
phenomenon  is  offered;  however,  a more  convincing  demonstration  would 
be  provided  by  using  monopropellaiit  componerts  differing  more  in  density 
than  AP  (density  = 1.95)  and  RDX  (density  =1.8).  It  would  also  be 
interesting  to  determine  for  mixture  containing  two  monopropcllants 
and  a fuel. 


TABLE  5.  Threshold  Conditiiins  for  RDX/AP  Mixtures. 


% AP 

Porosity 

Pore  size, 

urn 

P* 

MPa 

g/ cm2-s 

rag/ cm-s 

Wt. 

Vol. 

100 

0.  465 

119 

1.20 

0.315 

3.74 

75 

73.5 

0.  457 

115 

0.  35 

0.  350 

4.05 

50 

48 

0.450 

L15 

0.  28 

0.  553 

6.35 

25 

23.5 

0.  438 

107 

0.  50 

0.607 

6.  48 

0 

0 

0.426 

106 

1 . 20 

0.630 

6.65 

In  the  tests  on  various  oxidizer/fuel  combinations,  it  was  found 
that  the  threshold  conditions  for  onset  of  convective  burning  (1*)  were 
determined  primarily  by  the  oxidizer  with  little  distinction  between 
melting  and  nonmelting  fuels  and  negligible  effect  of  parLicU?  size  or 
fuel  fractions  up  to  55%.  Tlie  value  of  for  compositions  containing 
the  monopropellant  AP  varied  only  sligljtly  from  the  value  of  3.7  mg/cm-s 
obtained  for  pure  AP . This  near  con.stant  value  of  <p^  for  AP  was  obtained 
even  though  the  flame  temperature  of  the  AP/sucrose  mixtures  varied  with 
composition  by  a factor  of  two.  The  higher  value  of  ==  19  nig/cm-s  for 
the  KP  mixtures  agrees  well  with  the  stability  limit  (highest  value  of 
(()  for  which  conductive  burning  occurs)  for  secondary  cxplt>sives.  'i'lu' 
statement  was  made  that  no  explanation  could  be  found  for  the  flifferent 
values  of  for  AP  and  KP . It  was  concluded  that  fuel  melting  is  not 
important,  liowever.  It  does  not  necessarily  follow  that  oxidizer  melting 
has  no  effect  on  burning  stability. 

Taylor  (1962a)  also  introduces  the  melt  layer  concept  to  i-xplain 
the  combustion  stability  of  fusible  propellants.  He  points  out  that,  at 
constant  bomb  prtussure,  the  material  with  the  lower  melting  point  should 
have,  a thicker  melt  Layer,  hence  s'uould  develop  convective  hiirning,  at  a 
higher  pressure.  He  cites  HMX  (m.p.  = 278'’C)  which,  has  n tliri'shold 
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pressure  of  1.3  MPa  compared  Lo  PEITI  (ni.  p.  = I40“C)  with  a Llireshold 
pressure  of  2.8  MPa.  Andreev  (1966)  refutes  Taylor  with  three  c’ountcr 
arguments : 

1.  Andreev  believes  that  the  melt  layer  thickness  does  not  depend 
simply  on  melting  but  upon  boiling  point  as  well  (UMX  b.p.  > 340°C;  PETN 
b.p.  = 270°C) . 

2.  Taylor  explains  the  stability  <.)!  conductive  burning  of  NC  by  a 
melt  layer.  Andreev  dues  not  agree  that  NC  melts. 

3.  Andreev  claims  that  at  constant  pressure,  loss  of  stability  of 
conductive  bi  ning  at  iilgh  pressures  is  a restil  L of  the  higlier  dynamic 
pressure  which  is  generated  at  higher  static  pressures.  Taylor'. s re- 
joinder is  that  oombu.stion  of  secondary  explo.sives  cannot  generate  a 
sufficiently  high  dynamic  pressure  to  justify  Andreev's  argument. 

It  will  be  seen  (Section  3. 2. 2. 2)  that  tlu'  question  of  tin.'  role  of 
the  melt  layer  under  conditions  of  rising  pressure  (Ul.DCB)  i.s  no  more 
resolved  than  for  the  constant  pressure  condition  just  reviewed. 

Thus  far  in  this  section  only  criteria  neces.sary  for  onset  of 
convective  burning  have  been  ox.amined  along  with  the  effect  o1  particle 
size  and  a limited  range  of  compositions.  Even  less  information  is 
available  concerning  tite  rate  of  propagation  of  convective  burning  into 
porous  beds. 

Experimental  observations  reveal  that  convective  hu...ing  under 
constant  pressure  (spontaneous)  proceeds  relatively  smoethly,  ch.aracter- 
Ized  by  a plane  combustion  front  with  ,a  definable  velocity.  On  the 
other  hand,  convective  combustion  driven  by  a substantial  pressure 
difference  (forced  penetration),  such  as  i.s  found  in  tlu'  manomelric  bomb 
or  in  the  "embedded"  charge,  is  marked  by  irrep.ularities  resulting  in 
tongues  of  reaction  extending  deep  into  the  porous  m.aLerials.  The 
remainder  of  this  section  will  present  experimental  results  relevant  to 
the  rate  of  penetration  of  eom'oustlon  under  spontaneous  conditions. 

Figure  26  [Bobolev  (1966)],  used  earlier  to  illustrate  the  effect 
of  particle  size  on  tItreslioLd  pressure,  also  shows  the  depcndencs'  of 
the  resulting  convective  mass  burning  rate  of  RDX  on  pressure  and 
particle  size.  Since  a ne.arly  consL.ant  charge  density'  was  used  in 
obtaining  the  data,  tlie  stci-p  convective  burning  rate  curves  are  also 
representative  of  the  linear  velocity  ot  propa;;ation  of  Llie  combustion 
front  through  tluf  test  sample.  Although  the  slopes  of  the  convective 
burning  rate  curves  appear  to  increasi'  regularly  with  decreasing  par- 
ticle size,  scatter  and  paucity  of  data  ohviato  any  g.eneral  conclusion. 
The  important  feature  reve.aled  by  this  data  is  the  gre.ntlv  increased 
.sensitivity  of  m.ass  bu'^ning  rate  to  ch.anges  in  p ri'ssii ri‘  coiiipari-d  to  the 
conductive  mode. 
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Andreev  and  Chuiko  (1963)  studied  the  convective  burning  of  PETN, 
tetryl  and  NC  at  pressures  up  to  100  MPa  and  for  several  porosities  and 
particle  sizes.  Results,  shovm  in  Figure  28,  are  consistent  with  those 
for  RDX;  namely,  stability  of  conductive  burning  increases  with  decreas- 
ing particle  size  and  decreasing  porosity.  For  large  particles  (Curves 
1 and  2),  the  decrease  in  slope  at  high  pressures  is  explained  by  the 
increased  mass  of  inert  pore  gases,  hence  greater  dilution  at  higher 
pressures.  The  more  complicated  nature  of  Curves  3,  4,  and  5 (small 
particles)  is  attributed  to  competing  effects  of  four  phenomena: 

1.  Pore  gas  dilution  effect  is  greater  at  high  porosities  and 
pressures. 

2.  Permeability  to  gas  flow  is  greater  at  high  porosities  and  for 
large  particles. 

3.  Combustion  is  more  complete,  nearer  the  surface  for  high 

pressures,  hence  the  penetrating  gases  represent  a more  ener- 
getic ignition.  i 

4.  The  , dynamic  (combustion  generated)  pressure  is  greater  at  high 
burning  rates  and  pressures. 


For  the  high  porosity  (72%)  represented  by  Curve  3,  the  dilution 
effect  is  noticeable  even  at  comparatively  low  pressures;  for  the  small 
particles  (5  ym),  dilution  is  not  overcome  even  by  the  high  pressure 
effect  of  increasing  flame  temperature  and  dynamic  pressure.  At  lower 
porosity  (Curve  4)  the  effect  of  dilution  does  not  appear  until  a higher 
pressure  is  reached.  At  still  higher  pressures,  the  dilution  effect  is 
overcome  by  the  higher  combustion  temperature  and  dynamic  pressure.  At 
the  lowest  porosity,  the  dilution  effect  is  not  obvious  over  the  range 
of  pressures  studied.  This  competition  among  several  effects  is  illus- 
trated more  graphically  in  Figure  29,  which  shows  the  mass  burning  rate 
of  PETN  as  a function  of  relative  density,  6,  for  5 ym  particles  and  100 
MPa  pressure.  The  maximum  in  the  curve  arises  from  the  tradeoff  between 
dilution  and  permeability.  The  rising  portion  of  the  curve  is  attribut- 
able to  the  decrease  in  dilution  as  6 increases  while  the  declining 
portion  is  caused  by  the  decrease  in  permeability  and  penetration  rate 
at  high  densities.  Although  further  information  is  not  available,  it  is 
expected  that  relationships  similar  to  Figure  29,  would  be  obtained  for 
different  pressures,  particle  sizes,  and  materials,  depending  on  the 
interplay  among  the  four  effects  mentioned  above.  It  would  be  interest- 
ing to  investigate  (1)  whether  a maximum  mass  burning  rate  always  exists, 
as  in  Figure  29,  and  (2)  the  character  of  the  relation  at  rocket  motor 
pressures  for  multimodal  particle  sizes  and  standard  propellant  formula- 
tions. .A  p.innif leant  result  for  the  tests  deseribed  by  Andreev  end 
Chuiko  (1967,)  is  that  transition  to  detonation  never  oeeurred  eveti  at 
pressures  up  to  100  IlPa,  although  hurninj  rates  exceeded,,  coyuiuctivc 
rates  h:j  fistors  rancinn  from  10  to  120.  A possible  explanation  is  tliat 
the  size  of  the  test  sample  (5  mm  diameter  by  35  mm  length)  was  outside 
the  range  of  threshold  diameter  and  runup  length  for  the  materials 
studied. 
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FIGURL  28.  Pres.suro  Dependence'  of  Mass  Burning 
Rate  of  PHTN  at  Different  Relative  Density  and 
Particle  Size.  [Data  of  Andreev  and  Chuiko 
(1963),  taken  from  Belyaev  (1973),  Fip/are  60.] 
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FIGURE  29.  Dependence  of  Convective  Mass  Burning 
Rate  on  Relative  Density  for  PETN  (100  MPa,  Parti- 
cle Radius  20  um) . [Data  of  Andreev  and  Chuiko 
(1963)  taken  from  Belyaev  (1973),  Figure  61.] 


Taylor  (1962a)  conducted  a series  of  tests  to  determine  the  end 
effect  of  the  porous  bed  on  propagation  of  the  convective  burning  front 
in  HMX.  Conditions  which  remained  the  same  throughout  were:  particle 
size  range,  200-600  pm;  porosity,  43.2%;  bed  diameter,  6 mm;  bed  length, 

72  mm;  and  pressure  2.7  MPa.  Progress  of  combustion  was  determined  by 
cinephotography.  The  relation  between  time  and  position  of  the  flame 
front  is  shown  in  Figure  30  for  three  different  end  conditions.  When 
the  end  is  closed  (Curve  c) , the  propagation  velocity  accelerates  from 
the  beginning  until  the  counter  pressure  produced  by  the  closed  end 
causes  a deceleration.  For  the  open  end  (Curve  b),  acceleration  proceeds 
unhindered  to  the  end  of  the  channel.  Curve  (a)  shows  the  increase  in 
velocity  provided  by  a flow  restriction  (nozzle)  and  consequent  pressure 
increase  at  the  burning  end  of  the  charge.  It  is  noteworthy  that  for 
the  constant  pressure  conditions  of  the  tests,  there  is  no  initial 
deceleration  of  inflowing  gases  as  was  reported  in  Section  3.2.1  dealing 
with  an  ’’embedded"  charge.  Limitations  of  available  data  do  not  allow 
an  e.xplanation  of  this  distinction  at  the  present  time. 

Btjiyaev  (1973)  presented  an  interesting  empirical  relation  between 
propagation  rate  of  convective  combustion  and  conditions  for  its  onset 
as  shown  in  Figure  31.  The  abscissa  (R,},)  is  the  ratio  where  b is 

the  actual  value  and  is  its  onset  value  (see  earlier  discussion  in 

this  section).  The  ordinate  (Ry)  is  the  velocity  of  the  convective 
front  divided  by  the  normal  conductive  burning  rate.  It  is  not  clear 
whether  mass  or  linear  velocities  are  involved  in  this  latter  ratio; 
porosities  are  nearly  constant  so  that  the  data  are  relatively  cons  i s L«.'n  t . 
It  is  si'cn  t!iat  R;  is  a measure  of  how  far  conditions  exceed  the  onset 
conditions  while  R^  measures  the  effect  on  propagation  rate.  Tliere  is 
som--  ; att-  r in  the  data;  liowever,  considering  tile  range  of  materials 
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FIGURE  30.  Distance  Time  History  of 
Convective  Burning  of  HMX.  [Data 
taker,  from  T.iylor  (1962a),  Figures  2 
and  3.  ) 


tHMX,  RDX,  FETN)  and  particle  .sizes  (20-400  ’.m),  the  correlation  is 
reasonable.  An  accompanying  theoretical  presentation  is  not  clear. 
Moreover  the  formulas  derived  do  not  explain  the  slope  break  indicated 
nor  do  they  reduce  to  R =1  when  R =1. 

V <}■ 

Ihis  concludes  tlu’  review  ot  paj'ers  on  convective  burning  in  blind 
l^orous  bed.s  under  conditions  of  con.stant  bomb  pres.sure. 
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FIGURE  31.  Correlation  of  Convective  Burning  Rate 
of  Propellants  With  Stability  Criterion.  [Data  taken 
from  Belyaev  (1973),  Figure  63.] 


3. 2.2. 2 Blind  Porous  Bed  Burned  in  the  HLDCB 

Requirements  for  onset  of  convective  burning  in  blind  porous  beds 
have  been  reported  in  five  references  [Belyaev  (1973),  Belyaev  (1966), 
Gorbunov  and  Andreev  (1967),  Bobolev  (1965b),  and  Andreev  and  Gorbunov 
(1963)]  covering  a wide  range  of  effects.  These  include:  propellant 
type,  porosity,  particle  and  pore  size,  permeability,  charge  diameter 
and  length,  and  melting  point.  Inasmuch  as  theory  and  experiment  both 
indicate  that  threshold  pressure  is  affected  by  charge  length  and  diam- 
eter, these  quantities  were  chosen  in  the  range  where  their  effects  were 
unimportant.  The  dependence  of  threshold  pressure  on  charge  diameter 
for  PETN  is  shown  in  Figure  32.  Similar  considerations  of  length  led  to 
the  choice  of  a specimen  having  a diameter  of  10  mm  and  a length  of  50- 
100  mm. 

A series  of  tests  was  conducted  [Belyaev  (1966)]  to  compare  mixed 
(composite)  and  secondary  (single  ingredient)  explosives — a bomb  loading 
density  of  0.05  g/cm^  was  used,  producing  a dp/dt  of  0.1-10  bar/ms.  The 
results  depicted  in  Figures  33  and  34,  show  the  dependence  of  threshold 
pressure  on  porosity  for  a series  of  materials  having  a particle  size 
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FIGURE  32.  Effoi't  ol  Clii)rj;o  Dianiotor  on 
Threshold  Pressure  for  PETN.  |DaCa  taken 
from  Belyaev  (1966),  Figure  8.] 
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FIGURE  33.  Threshold  Pressure  as  a 
Function  of  Porosity  for  Secondary 
Explosives  (Particle  Radius,  10-20 
um) . [Data  taken  from  Belyaev  (1966), 
Figure  3. ] 
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(1)  Bltumen/KCl04  (3)  TNT/AP 

(2)  Polystyrene/KClO^  (4)  Bitumen/AP 

FIGURE  34.  Threshold  Pressure  as  a 
Function  of  Porosity  for  Mixed  Ex- 
plosives (Particle  Radius,  10-20  ym) . 
[Data  taken  from  Belyaev  (1966) , 
Figure  14. ] 


range  of  10-20  ym.  It  is  noteworthy  that  threshold  pressures  are  much 
lower  for  the  composite  materials^  indicating  their  greater  susceptibil- 
ity to  convective  burning.  Belyaev  interprets  the  results  shown  in 
Figure  34  as  indicating  that  development  of  convective  burning  is 
governed  primarily  by  the  oxidizer  properties.  While  there  may  be 
unreported  evidence  to  support  this  claim,  it  should  be  noted  that  the 
results  presented  in  Figure  34  suggest  similar  effects  owing  to  changes 
in  oxidizer  or_  fuel.  The  important  conclusion  of  Figures  33  and  34  are 
that  for  otherwise  equivalent  situations:  (1)  the  threshold  pressure 
decreases  with  increased  porosity,  becoming  more  dependent  at  porosities 
below  0.05;  (2)  the  threshold  pressure  is  highly  sensitive  to  propellant 
type  at  low  porosity  and  much  less  dependent  at  porosities  exceeding 
0.15.  The  implication  is  that  at  the  low  porosities  most  often  en- 
countered in  practice,  the  propellant  properties  have  a large  influence 
on  development  of  convective  burning.  Stated  somewhat  differently; 
Figure  shows  that  homogeneous  (single  component)  systems  have  a 
relatively  high  threshold  pressure  while  Figure  34  shows  that  composite 
(heterogeneous)  propellants  have  pressures  that  are  of  concern  when 
compared  to  rocket  ahofnber  pressures  where  the  porosities  lie  also  in 
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the.  area  of  anneer'n.  At  higher  porosities  (>0.15)  this  effect  decreases 
until  propellant  properties  have  only  a minor  effect.  The  curves  of 
Figures  33  and  34  are  approximated  by  the  empirical  relationship: 

p^(m  - B)  = A (11) 

where 

p^  = threshold  pressure 
m = porosity 

A,  B = constants  dependent  on  propellantj  with  B representative  of 
unconnected  porosity 

The  foregoing  results  are  for  propellants  with  a particle  size 
range  of  10-20  pm.  The  results  of  a series  of  tests  on  PETN  for  dif- 
ferent initial  particle  sizes  are  plotted  in  Figure  35.  It  is  seen  that 
large  particles  contributed  to  a lower  threshold  pressure  and  that  the 
effect  of  particle  size  is  greater  at  low  porosity. 


Particle  radius,  urn 

(1)  10-20 
(2)  120 
(3)  550 


FIGURE  35.  Effect  of  Porosity  and 
Particle  Radius  on  Threshold  Pres- 
sure for  PETN  (Charge  Diameter, 

5 mm).  [Data  taken  from  Belyaev 
(1966) , Figure  5. ] 


I 

i 


A 


A 

1 

.1 


51 


NWC  TP  6007 


The  dependence  of  threshold  pressure  on  gas  permeability  was  also 
investigated.  The  permeability,  k,  is  defined  by  the  equation: 

-1/A^  dU/dt  ■ (k/u)(8p/3x)  (12) 

where 


Af,  =■  sample  cross-sectional  area  normal  to  gas  flow 
U = volumetric  throughput 
k =*  sample  gas  permeability 

and  represents  the  volume  flow  rate  per  unit  area  for  a unit  fluid 
viscosity  and  pressure  gradient.  For  A^.  in  cm^,  dU/dt  in  cm^/s,  y in 
centipoise  (cP),  and  8p/3x  in  atm/ cm,  k is  given  in  darcies.  Figure  36 
shows  the  relationship  between  permeability  and  porosity  for  several 
propellants  and  particle  sizes.  The  influence  of  propellant  type  is  not 
explained  but  is  probably  a result  of  different  particle  size  distribu- 
tion for  the  same  range  of  particle  sizes.  If  selected  data  from  Fig- 
ures 33  through  35  are  cross-plotted  with  data  from  Figure  36,  the 
effect  of  permeability  on  threshold  pressure  is  obtained  as  shown  in 
Figure  37.  The  composite  propellants  are  well  correlated  by  the  equa- 
tion 


k exp(p^/D)  = C (13) 

where  C and  D depend  weakly  on  the  propellant  formulation.  The  effect 
of  particle  size,  other  conditions  remaining  the  same,  is  shown  in 
Figure  38  for  PETN.  The  importance  of  Figures  3"  and  38  is  that: 

(1)  secondary  explosives  such  as  TNT  and  PETN  are  less  susceptible  to 
convective  burning  and  have  a more  complicated  and  sensitive  dependence 
of  threshold  pressure  on  permeability  than  do  the  composite  propellants; 

(2)  the  dependence  of  threshold  pressure  on  propellant  properties  in- 
creases with  decreasing  porosity;  and  (3)  permeability j like  porosity 
and  particle  size\  does  not  uniquely  determine  critical  pressure.  An 
explanation  offered  for  the  failure  of  permeability  to  account  for 
breakdown  of  surface  burning  is  that  burning  penetrates  the  larger  pores 
first  but  that  permeability,  a parameter  which  integrates  the  flow 
properties  of  all  sizes  of  pores,  does  not  reflect  pore  size  distribu- 
tion. An  important  deficiency  in  these  studies  is  the  failure  to  estab- 
lish the  effect  of  the  rate  of  pressure  rise  on  threshold  pressure. 

While  the  results  just  described  were  obtained  with  a geometrically 
blind  porous  bed,  care  was  taken  to  choose  a bed  length  so  that  the 
presence  of  the  closed  end  did  not  affect  the  results.  A difference 
between  geometric  restriction  versus  hydrodynamically  blind  needs  to  be 
kept  in  mind.  That  is,  the  test  bed  was  made  long  enough  to  avoid 
geometric  restriction.  In  additional  auxilliary  tests,  the  end  effect 
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0 O.Ob  0.10  0.10 


POROSIlv,  ,n 

Particle  radius,  pm 

(1)  PETN  (10-20)  (3)  TNT  (10-20) 

(2)  TNT/AP,  stoichiometric  (10-20)  (.4)  PKTN  (.3601 

FIGURE  36.  Effect  of  Porosity,  Particle  Size, 
and  Material  on  Permeability.  [Data  taken  from 
Belyaev  (1966),  F.igure  6.] 
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(1)  TNT  (4)  TNT/AP 

(2)  PETN  (3)  Bitumen/AP 

(3)  Bl  tunien /KCl O, 

A 

FIGURE  37.  Dependence  of  Threshold  Pressure  on 
Permeability  and  Material  for  Particle  Radius 
of  10-20  im;.  [Data  taken  from  lielyaev  (1966), 
Figure  7a. ] 
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(1)  10-20  urn 

(2)  550  ym 

FIGURE  38.  Threshold  Pressure  for  PETN 
as  a Function  of  Permeability  for  Two 
Particle  Sizes.  [Data  taken  from  Belyaev 
(1966),  Figure  7b.] 


was  found  to  be  negligible  if  (1)  the  permeability  was  less  than  10 
darcy,  and  (2)  at  least  15-30  mm  of  unburned  length  remained  at  the  time 
was  attained.  For  high  permeability  and  short  beds,  the  threshold 
pressure  is  stated  to  be  sensitive  to  length  [Belyaev  (1966)],  but  no 
quantitative  results  are  presented. 

Andreev  and  Gorbunov  (1963)  determined  the  threshold  porosity  at 
which  convective  burning  developed  for  TNT,  picric  acid,  PETN,  RDX,  and 
mercury  fulminate. ^ These  results,  along  with  flow  resistance  (propor- 
tional to  reciprocal  of  permeability)  are  indicated  in  Table  6.  The 
extent  to  which  the  quantities  of  Table  6 are  comparable  is  uncertain, 
since  charge  weights  and  dimensions  varied  among  materials.  An  HLDCB 
was  used  having  a volume  of  50  cm3  with  an  igniter  consisting  of  1 g of 
black  powder  which  produced  an  initial  pressure  of  50  atm  (5  MPa) . 
Despite  the  poorer  control  of  conditions,  results  show  a ranking  of 
susceptibility  to  convective  burning  similar  to  that  indicated  in  Figure 
33.  The  greater  stability  of  PETN  over  RDX  indicated  in  Table  6 may  be 
due  to  the  smaller  particle  size  of  the  PETN.  Additional  experiments 
were  conducted  on  PETN  to  show  the  effect  of  charge  length  and  porosity. 
The  change  from  normal  to  convective  burning  as  charge  length  increased 
(Figure  39)  indicates  the  stabilizing  effect  of  the  closed  or  blind  end 
and  is  in  quaUtativc  agreement  with  theory  (see  section  on  theory). 
Figure  'lO  shows  that  convective  burning  is  more  likely  to  occur  early  in 
tli>'  u-st  (al  a lower  pressure)  if  a large  particle  size  is  used.  An 
Lr'j'or*  iru  f'-ature  described  for  the  tests  whose  pressure  time  records 
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TABLE  6.  Comparative  Stabilities  of  Conductive 
Burning  of  Explosives 


1 j 

Material 

Threshold 

porosity 

Threshold  flow 
resistance  (relative) 

TNT 

0.  33 

2.5 

Picric  acid 

0.19 

21 

RDX 

0.07 

110 

PETN 

0.  18 

400 

Mercury  fulminate 

0.  11 

1900 

Par  t ic  1l‘ 
size,  um 

50-60 

10-20 

50-60 

10 

50-60 
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Particle  size 

(1)  3O0-A00  and  460-500  utn 

(2)  100-160  um 

(3)  63-100  ym  (detonated) 

(4)  10  ym 

FIGURE  40.  Effect  of  Particle  Size 
on  Burning  History  of  PETN  (Constant 
Length  and  Porosity).  (Data  taken 
from  Andreev  and  Gorbunov  (1963), 
Figure  7 . ] 


are  shown  in  Figure  40  is  that  the  most  violent  combustion  occurred  for 
particles  of  Intermediate  size  where  an  optimum  trade-off  exists  between 
permeability  (larger  for  large  particles)  and  specific  surface  area 
(larger  for  smaller  particles). 

One  explanation  for  the  greater  stability  of  conductive  burning  of 
secondary  explosives,  e.g.,  TNT,  picric  acid,  DINA,  PETN,  and  RDX  (in 
order  of  decreasing  stability)  is  that  these  materials  melt  during 
burning,  forming  a liquid  surface  layer  which  blocks  the  flow  of  gases 
into  defects  [Belyaev  (1966)].  The  criterion  advanced  by  Belyaev  for 
onset  of  convective  burning  is  the  discontinuity  of  the  surface  liquid 
layer,  which,  in  turn,  requires  that  the  melt  layer  thickness  be  less 
than  the  maximum  pore  size.  The  melt  layer  hypothesis  is  supported  by 
the  fact  that  calculated  melt  layer  thicknesses  of  the  five  materials 
listed  above  give  values  which  decrease  in  the  same  order  as  the  ob- 
served stability.  An  attempt  was  made  to  explore  this  phenomenon 
[Gorbunov  and  Andreev  (1967)]  by  testing  pairs  of  chemically  similar 
materials  with  significantly  different  melting  points  for  their  propen- 
sity to  develop  convective  burning.  Particle  sizes  were  held  constant 
at  jO-100  cm.  Results,  shown  in  Table  7,  are  to  d Mnoiis r r.a 

that  the  melt  layer  does  not  aihord  protection  against  onset  of  convec- 
tive burning. 
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TABLE  7.  Effect  of  Melting  Properties 
on  Stability  of  Conductive  Burning. 


Relative  threshold 
resistance 

Threshold 
porosity,  % 

Me  1 1 in  g 
point,  °C 

TNT 

2.8 

33 

80 

Xylyl 

3.0 

27 

178 

Trinitro  benzene 

29 

12 

121 

Hezanitro  biphenyl 

26 

IR 

238 

RDX 

no 

7 

203 

HMX 

180 

6 

272 

DINA  (15°C) 

55 

13 

49 

DINA  (45°C) 

55 

12 

49 

In  addition  to  tho  effect  of  melting  point,  the  effect  of  initial 
temperature,  hence  melt  layer  thickness,  was  investigated  for  DINA. 
Results,  also  shown  in  Table  7,  suggest  that  the  thickness  of  the  imper- 
meable liquid  layer  does  not  govern  onset  of  convective  burning.  Numer- 
ous phenomena  such  as  reactivity  of  tho  liquid  layer,  boiling  point, 
surface  tension,  and  viscosity  were  not  investigated.  Also,  results 
discussed  in  this  section  of  the  report  are  for  rapidly  rising  pressures 
Tlie  effect  of  a melt  i;iycr  in  a near  constant  pressure  environment  was 
reported  in  Section  3. 2. 2.1. 

The  final  paper  dealing  with  porous  beds  under  conditions  of  rising 
pressure  [Bobolev  (1965b)],  considers  the  effect  of  addition  of  paraffin 
to  RDX  (phlegmatl zatlou) . Comparison  was  made  of  burning  of  pure  RDX 
(particle  size,  200  urn)  and  phlegmatizcd  RDX  (particle  size  range,  50- 
360  urn).  Results  are  shown  in  Figure  41  for  levels  of  paraffin  varying 
from  zero  (pure  RDX)  to  107,.  For  this  series  of  experiments,  the  addi- 
tion of  paraffin  has  no  significant  effect  on  tho  relation  between 
relative  density  and  breakdown  pressure;  however,  even  0.5%  paraffin 
prevents  the  buildup  of  convective  burning  Into  detonation.  The  essence 
of  the  explanation  is  that  the  paraffin  acts  botfi  as  a heat  barrier, 
preventing  heat  transfer  to  the  active  com|)onent  (RDX),  and  as  a cooluit 
(endothermic  melting)  for  the  combustion  products.  Thus,  the  role  of 
the  phlegmatizer  is  confined  to  stages  of  DDT  following  onset  of  convec- 
tive burning,  involving  accelerated  combustion  and  formation  of  weak 
shocks . 

This  oiiv.  ludes  thic  tevi  cw  of  c xp c — inr— n t a i p a c r s on  a u r . > 1 1 . g in 
blind  |)orous  beds  under  conditions  of  risii\g  pressure. 
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(1)  Pure  RDX 

(2)  10%  Paraffin 

FIGURE  41.  Effect  of  Addition 
of  Paraffin  on  Threshold  Pres- 
sure of  RDX  (Particle  Radius, 
200  urn).  [Data  of  Bobolcv 
(19656)  taken  from  Belyaev 
(1973) , Figure  32.  J 


3 . 2 . 3 . 1 Open  Porous  Bed  Burned  in  the  LLDCB 

Taylor  (1962a)  determined  the  effect  of  pressure  and  partlc^li'  size 
on  mass  burning  rate  of  HMX . His  samples,  which  were  burned  in  n 
Crawford  bomb,  consisted  of  6 mm  diameter  paper  tubes  of  7.8  cm  length 
filled  with  the  granular  material.  A plug  of  plastic,  cemented  into  the 
lower  end  of  the  tube  to  retain  the  sample,  wan  not  believed  to  provide 
a gas  tight  seal.  Experimental  results  are  shown  in  Figure  42  which 
shows  no  convective  burning  for  small  particles  (5  nm)  even  at  liigh 
pressures.  Although  bulk  densi..y  (porosity)  varied  somewhat  (1.02-1.20 
g/cm3)  throughout  the  particle  size  range,  it  is  clear  that  (1)  onset 
of  convective  burning  shifts  to  lower  pressures  for  large  panicles,  and 
(2)  convective  burning  rate.s  are  at  Jc.ist  an  order  of  magnitude  grtater 
than  conductive  and  much  more  sensitive  to  pressure.  Similar  trends 
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FK'iURF,  42.  Convective  Mass  Burning  Rati.' 
of  HMX  ns  a Function  of  Particle  Size  and 
Pre.ssure.  [Data  taken  from  Taylor  (1962a), 

Figure  1 . ] 

have  been  noted  and  reported  for  otlier  experimental  conditions  and 
materials.  Effect  of  the  end  condition  has  been  presented  in  Siction 
3. 2. 2.1  (see  Figure  30)  and  shows  the  effect  of  the  closed  end  in  making 
development  of  convective  burning  more  difficult. 

In  a separate  Investigation,  Taylor  (1962b),  using  similar  tecli- 
nirjues,  measured  the  effect  of  pressure  and  HMX  particle  sizi'  on  mass 
burning  rate,  Tiie  effect  of  pressure  i.s  shown  in  Figure  43  for  two 
different  combinations  of  density  and  jiarticle  size.  Not  only  is  mas.s 
burning  rate  nearly  independent  of  particle  size  (density  varying)  but 
conductive  burning  prev.’iils  ovi-r  the  entire  pre.ssure  range,  as  may  be 
verified  by  comparison  with  the  6 nm  curve  of  Figure  42.  At  the  same 
time,  the  linear  burning  r^Lt*  v.irie.s,  approximate  Iv  inversely  witii 
density  (density  varying  from  1.02-1.66  g/iin^).  Tin.'  el  feet  ot  particle 
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FIGURE  43.  Efft-ct  of  Pressure,  Density, 
and  Particle  Size  on  Mass  Burning  Rate  of 
HMX  (Conductive  Region).  [Data  taken  from 
Taylor  (1962b),  Table  2.) 


size  on  mass  burning  rate  is  sliown  in  Figure  44  for  the  relatiVL‘ly  low 
pressure  of  12.6  atm  (1.26  MPa)  where  convective  burning  doi-s  not  di‘v- 
elop  over  the  range  of  conditions  studied.  It  is  seen  that,  aitliough 
convective  burning  is  absent,  the  mass  burning  rate  increases  sliglitly 
(compared  to  effect  of  particle  size  on  convective  burning)  with  ]iarti- 
cle  size.  Similar  results  were  obtained  with  PETN,  showing  that  in  tlu' 
conductive  mode,  mass  burning  rate  is  relatively  insensitive  to  particle 
size  and  density. 

Two  alternate  explanations  are  offered  for  the  above  phenomena. 
Taylor  attributes  the  insensitivity  of  mass  burning  rate  to  particle 
size  to  a smooth  melc  layer  which  comprises  a pressure  dependent  source 
of  material  fed  to  the  gas  phase.  It  is  only  at  pressures,  hence  burn- 
ing rates  and  melt  layer  thicknesses,  which  permit  the  melt  layiT  tc> 
conform  to  surface  irregularities  without  losing  continuity  that  the 
burning  rate  is  affected  by  pore  size.  At  even  higher  pre.ssures,  tin' 
protection  afforded  by  the  melt  layer  disappear.s,  leading  to  convective 
burning . 
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FIGURE  44.  Kflect  ol'  P.irticU'  Size 
on  Mass  Burning  Rato  of  UMX  at  1.26 
MPa  (Conduetive  Region)  for  Nearly 
Const.int  Density.  (Data  taken  from 
Taylor  (1962b),  Table  4.] 

Andreev's  explanation  is  more  tenuous  and  less  complete,  depending 
an  his  claim  that  heat  transfer  from  the  gas  to  the  solid  surface  is 
independent  of  porosity.  Without  invoking  a melt  layer,  Andreev  ex- 
plains the  results  on  the  basis  of  gas  penetration  which  would  occur  .it 
a higher  rate  for  larger  particles.  The  issue  is  not  yet  resolved. 


3.3  DAMAGED  PROPELLANT  BURNED  IN  THE  LLDCB 

The  cornbustion  of  damaged  propellants  was  studied  at  the  Nava] 
Weapons  Center.  Two  different  series  of  tests  were  performed:  one  in 
whicli  the  propellant  was  strained  before  ,and  during  tlie  combustion  and 
one  in  which  the  propellant  was  strained  until  severe,  damagi'  oeeurred 
but  the  imposed  force  wa.s  removed  prior  to  the  sample's  bi'ing  burned. 
Boggs,  Zurn,  and  Derr  (1976)  studied  the  first  ease.  Their  results 
(Figures  45  and  46)  showed  some  augmentation  of  burn  rate  at  some  values 
of  strain  and  pressure. 

Two  propellants  were’  used:  a composite  propellant  with  inert 
rubber  ba.se  hinder  (Figure-  45),  and  a high  enorj-y  cross-linked  doubi-  - 
base  propellant  (Figure  46).  These  propellants  were  placc-d  in  sample 
holders  as  shoc-m  in  Figure-  47.  A .small  piece-  of  microscope  slide  was 
plaee-d  against  tlu^  sample  sideji  to  pre-vent  the’  cold  bomb  gases  ironi 
filling  the  voids  produee-d  by  the  subseejut-nt  .straining  accump  1 i shetl 
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FIGURE  47.  Sample  Holder  Used  in 
Measuring  Conductive  Burning  Rate 
of  Strained  Propellants. 


using  the  device  of  Figure  47.  The  sample,  under  strain  in  the  fixture, 
was  burned  in  a window  bomb  (loading  density  less  than  0.001  g/cm^). 

Bum  rates  were  determined  from  the  films. 

Data  for  both  types  of  propellants  (composite  and  cross-linked 
double-base  high  energy  propellants)  show  burn  rate  increase  at  pres- 
sures and  strains  greater  than  some  threshold  values.  It  was  emphasized 
that  both  threshold  values  have  to  be  exceeded,  exceeding  just  one  was 
not  sufficient.  For  example  high  strain  but  low  pressure  did  not  cause 
augmentation  nor  did  high  pressure  but  low  strain. 

The  mechanical  response  of  the  propellants  to  strain  was  studied 
using  a binocular  microscope.  These  studies  showed,  using  Propellant 
A as  an  example,  that  at  4%  strain,  debonds  (separation,  on  a micro-scale, 
of  the  solid  particle  from  the  polymeric  binder)  between  ingredients 
occurs.  Between  9-11%  strain,  these  debonds  are  often  fully  developed 
cracks,  with  the  walls  of  the  crack  in  close  proximity.  At  approximately 
16%  these  cracks  are  open  voids;  that  is,  the  walls  of  the  crack  are  no 
longer  in  contact  with  one  another.  At  approximately  24%  the  sample  is 
often  riddled  with  large  cracks  and  the  sample  fails. 

The  authors  provided  a mechanistic  explanation  for  the  burn  rate 
augmentation  due  to  strain  and  pressure.  At  low  strain  values  the 
propellant  was  not  significantly  damaged  and  so  regardless  of  the  flame 
stand-off,  augmentation  did  not  occur.  When  the  propellant  was  highly 
strained  and  fissured,  augmentation  occurred  if  the  flames  penetrated 
into  these  fissures.  At  low  pressures  the  flame  stood  too  far  from  the 
surface  to  allow  penetration,  but  at  high  pressures  the  flame  was  close 
enough  to  the  surface  to  penetrate  the  fissures  and  cause  burn  rate 
augmentation. 
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Tlic  abovo  ti'sL.s  wore  ooiiituotod  on  samples  with  imposed  strain  (rhal 
is,  the  voids  were  open).  Another  series  of  tests  was  conducted  cn  pro- 
pellant samples  which  were  heavily  damap.ed  but  the  imposed  strain  was 
relea.sed  prior  to  burning.  Samples  of  propellant  were  subjected  to 
tensile  strains  whie.h  caused  significant  damage.  The  sample  was 
strained  to  25%  and  then  allowed  to  relax  to  zero  .strain.  The  sample 
was  then  strained  to  approximately  35Z  where  the  st ress/strain  character- 
istics indicated  severe  damage.  Thesi’  d;.)maged  samples  woi'i-  burned  at 
pressures  up  to  34.5  MPa.  The  re.suits  are  presented  in  Figure  48.  It 
was  found  that  damaged  and  undamaged  propel  l.ints  h.ad  the  same  burn  rate, 
indicating  that  damaged  pnoellant  tends  to  burn  normally  if  tie  strain 
induced  voids  are  allowed  to  close  (a  debond  condition  existed,  but 
there  was  no  open  pore) . 


FIGURK  48.  Effect  of  Damage  on  Conductive 
Burning  Rate  of  a High  Energy  Cross-Linked 
Double-Base  Propellant. 


4.0  THEORY 

As  outlined  earlier  (Section  2.2),  the  literature  on  the  simplified 
theory  of  oarty  convective  burning  falls  Into  two  principal  categories: 
.single  pore  and  porous  bed.  Further  subdivision  is  according  to  the 
three  requirements  for  convective  burning,  viz.,  fluid  dynamic,  heat 
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transferj  and  propagation  of  an  ignition  front.  The  first  and  second 
subsections  which  follow  present  the  available  analyses  relevant  to 
single  pore  and  porous  bed  situations.  In  addition,  a third  section  is 
included  of  miscellaneous  results  which  do  not  fit  into  the  general 
morphological  pattern. 


4.1  SINGLE  PORE 

All  three  requirements  for  single  pore  convective  burning  have  been 
analyzed  theoretically  using  simplified  concepts.  The  fluid  dynamic 
criterion  considers  only  a blind  pore  with  a rising  external  pressure 
and  is  based  on  conservation  of  mass;  hence  frictional  effects  are 
ignored.  The  criteria  of  ignition  and  propagation  of  ignition  include 
a simplified  energy  balance  but  again  omit  the  momentum  equations. 

4.1.1  Fluid  Dynamics 


The  fluid  dynamic  requirement  of  gas  penetration  into  a single 
channel  can  be  represented  by  the  expression 

V > r (14) 

8 

where  Vg  is  the  rate  of  gas  flow  into  the  channel  and  r is  the  normal 
(conductive)  regression  rate.  The  condition  is  necessary  but  insuf- 
ficient for  the  propagation  of  an  ignition  front  into  the  pore.  Since  r 
is  a reasonably  well-known  function  of  environmental  pressure  (less 
well-known  as  a function  of  dp/dt  and  gas  velocity)  and  propellant 
temperature,  it  remains  to  calculate  Vg  for  the  appropriate  conditions. 
Ideally,  this  calculation  would  involve  application  of  transient  conser- 
vation equations  and  would  require  knowledge  of  transient  friction  and 
heat  transfer  in  the  entrance  section  of  the  channel.  The  only  known 
simplified  approaches  [Sobolev  (1965a)  and  Belyaev  (1973)]  ignore  these 
complications  and  consider  the  situation  depicted  in  Figure  49.  A 
closed  end  channel  of  length  L and  unit  cross-sectional  area  is  initial- 
ly filled  with  gas  at  temperature  Tq  and  pressure  p.  External  gas  at 
temperature  Tg  and  pressure  p,  equal  initially  to  pore  gas  pressure, 
flows  in  under  the  action  of  an  external  pressure  rise  dp/dt.  The 
following  assumptions  are  made  in  the  analysis: 

1.  ■ The  incoming  gases  mix  with  the  pore  gases,  the  mixture  tem- 

perature remaining  constant  at  the  initial  pore  temperature, 
Tgo- 

2.  The  pore  is  not  "too"  wide.  This  would  seem  to  imply  one 
dimensional  flow  (no  recirculation). 

3.  Relaxation  time  of  pore  pressure  is  less  than  relaxation  time 
of  external  pressure.  In  other  words  the  pore  process  is 
quasi-steady . 
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FIGURE  49.  Simp  1 if i<.>d  Modol  of  Gas  Flow 
into  Blind  Chaunol  Under  Rising  Pressure 
Conditions. 


4,  Tlie  product  Ldp/dt  is  "not  too  large."  This  is  related  to 

assumption  3 since  the  relaxation  time  of  the  pore  increases 
with  L while  the  relaxation  time  of  the  external  pressure 
increases  with  decreasing  dp/dt.  An  expression  for  v can  be 
derived  by  equating  the  quantity  of  incoming  gas  ^ 


th.  = V 0 = V p/bT 

Jn  g g g g 


(15) 


to  the 


increase  in  gas  in  the  pore 


m 

pore 


L p/bT 

go 


(16) 


during  the  Lime  dt.  The  specific  gas 
tion  (16)  depends  upon  assumption  1. 
required  re.sult; 


constant  is  denoted  by  b.  Equa- 

Equating  m.  to  m provides  the 

in  pore 


V =•  (L/p)(dp/dt)(T  /T  ) (17) 

g g go 

Recalling  Eq . (14),  one  notes  that  the  hydrodynamic  condition  for  inflow 
is  facilitated  (subject  to  the  restrictions  mentioned)  by:  (1)  deei)i.;r 
channels,  (2)  low  pressures,  (3)  high  rate  of  pressure  rise,  (4)  hot 
combustion  ))roducLs,  .and  (5)  cool  propellant.  The  effects  of  Low  pro- 
pellant temperature  and  low  pressure  are  twofold:  they  not  only  in- 
crease the  left-hand  side  but  decrease  the  right-hand  side  of  l^q . (14), 

thereby  improving  conditions  for  inflow.  It  should  be  stressed  that  E<]. 
(17)  was  derived  for  renditions  simulating  a blind  port?,  open  al  the 
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ignited  end  to  a rising  pressure.  Simplified  expressions  giving  Vg  in  a 
single  pore  for  other  conditions  have  not  been  found  in  the  current 
survey. 

4.1.2  Ignition  of  the  Wall  of  a Single  Pore 

Satisfaction  of  the  inequality  (Eq.  (14))  is  no  guarantee  that  the 
inflowing  hot  gases  will  heat  the  pore  wall  sufficiently  to  achieve  a 
condition  for  self-sustained  reaction,  usually  called  ignition.  An 
analysis  to  establish  conditions  for  ignition  of  pore  walls  by  hot  gas 
flow  is  concerned  with  two  broad  problem  areas:  (1)  definition  of 
ignition,  and  (2)  calculation  of  conditions  resulting  from  the  hot  gas 
flow.  The  simplest  ignition  criterion  is  an  ignition  temperature  which 
is  unique  to  the  material.  Such  a concept  ignores  the  variation  of 
ignition  temperature  with  experimental  conditions  such  as  heating  his- 
tory and  nature  of  the  surrounding  atmosphere.  Zeldovich  (1942)  pointed 
out  that  attainment  of  a critical  surface  temperature  alone  is  insuf- 
ficient to  assure  ignition,  but  that  the  entire  temperature  profile 
(variation  of  temperature  with  distance  from  the  heated  surface)  should 
resemble  that  under  steady-state  burning  conditions  for  the  substance 
being  ignited.  It  is  convenient  to  characterize  the  thermal  profile  by 
the  depth  of  penetration  of  the  heat  introduced  at  the  surface;  it  is 
the  dual  requirement  of  surface  temperature  and  thermal  profile  depth 
which  serves  as  ignition  criteria  in  the  analysis  by  Margolin  and  Chuiko 
(1965). 


The  second  problem  area  in  pore  wall  ignition  is  the  calculation  of 
surface  temperature  and  thermal  profile  depth  in  the  pore  wall  which 
results  from  ^nf lowing  hot  gases.  To  simplify  the  calculations, 

Margolin  and  Chuiko  (1965)  introduced  several  approximations.  Actually, 
the  inflowing  gases  cool  gradually,  simulataneously  heating  the  wall. 

It  is  assumed  that 


1.  The  gas,  at  constant  temperature,  Tg,  gives  up  its  heat  to  the 
pore  wall  at  constant  temperature,  Tw  This  neat  transfer  takes  place 
over  a cooling  length,  L^,  and  with  constant  rate  of  heat  transfer. 

2.  All  flow  variables  such  as  temperature,  velocity,  and  gas 
density  are  timewise  and  spacewise  constant  over  the  distance,  L^. 

3.  The  gas  flow  is  laminar  and  satisfies  the  inequality 


RePr 


> 0.1 


for  which  the  Nusselt  number  is  constant,  given  by 


(18) 


Nu=hd/A  =3.66 
P g 


(19) 
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where  x is  the  distance  from  the  pore  entrance,  h is  the  film  coeffici- 
ent for  heat  transfer,  dp  is  the  pore  hydraulic  diameter,  and  Ag  is  the 
thermal  conductivity  of  the  gas.  An  essential  step  in  the  establishment 
of  Ignition  conditions  is  the  calculation  of  the  length  of  the  cooling 
section  L^.  In  unit  time,  the  heat  lost  from  the  gas  to  the  wall  is 


q,  =■  (v  - r)  p (•  (T  - T )Ttd  jk 
loss  g g g g w'  p 


(20) 


This  same  heat  must  be  transferred  across  the  gas  film 

1. 


trans 


Tid  L h(T 

pc  g 


T ) 
w 


(21) 


Equating  the  two  expressions  permits  solving  for  L . 

L =•  ,:)  c d (v  - r)/4h  = (v  - r)d  ^/ANua  (22) 

eggpg  g p g 

where  the  '.hermal  diffusivity  of  the  gas,  a^,  is  defined  as 

n =“  A /(I  c (23) 

g g g g 

The  next  step  in  applying  Zeldovich's  criteria  consists  of  the 
evaluation  of  several  quantities: 

1.  The  time,  i,  available  for  heating  the  pore  wall 
The  depth,  i,  of  the  heated  layer  at  the  time,  r 


3. 

A. 

5. 


The  depth.,  f , of  the  steady-state  thermal  layer 
The  time,  i ',  required  to  establish  the  steady-state  thermal 
layer  ^ 

The  time,  required  to  heat  the  pore  wall  from  the  Initial 

temperature,  T^ , to  the  ignition  temperature,  T^. 

The  Lime  available  for  hcatlug  the  pore  wall  is  simply  the  time 


required  to  consume  the  heated  length  of  the  pore,  approximated  by  L , 
.-'t  the  conductive  rate  r Hence  ^ 


b /r 
c 


(2A) 


The  thermal  penetration  depth  Is  actually  infinite  since  the  tem- 
perature profile  at  any  Instant  is  asymptotic  to  the  initial  temper- 


ature, Tq.  However,  at  the  time  t,  the  temperature  at  depth  (tsO 


1/2 


has  dropped  by  90%  of  the  difference  Ty-Tf,;  this  distance  is  commonly 
used  as  the  depth  of  the  thermal  layer.  Therefore  we  may  write  approx- 
imately 


, ,1/2 
(u  1 ) 
s 

(25) 

. ,1/2 
(n  T ) 

(26) 
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The  time  required  to  establish  the  steady-state  thermal  layer  is 
again  infinite,  with  a leasonable  approximation  being  provided  by 


i=t/r  (<17) 
s s 

which  states  that  the  steady-state  profile  will  be  attained  by  the  time 
the  surface  has  regressed  a distance  equal  to  tlie  steady-state  penetra- 
tion distance. 


The  time  required  to  heat  the  pore  wall  to  the  ignition  temperature? 
(T^)  is  obtained  by  equating  the  heat  flowing  across  the  gas  film  to  the 
heat  absorbed  by  the  solid,  assuming  spacewise  constant  temperature  in 
the  solid. 


hi.  d (T  - T )i^  = 

r'  n CJ  W 7% 


S S C 


T )ld  (a  7 
o p s 


,1/2  ^ 


(28) 


Solving  for  t*  gives 


where  the  subscript  s refers  to  solid  phase.  The  quantity  a is  a geo- 
metric parameter  with  the  value  of  zero  for  rectangular  slits  and  unity 
for  circular  pores,  while  P is  given  by 

B = A (T.  - T )/A  (T  - T.)  (30) 

s * o g g * 

Zeldovich'.s  criteria  are  given  by 

t/i^  > i (31) 

P/v,  >lort/T>l  (32) 

s .s 

Equation  (31)  expresses  the  requirement  that  the  available  pore  heating 
time  is  long  enough  to  heat  tlie  pore  wall  to  the  ignition  temperature. 
Equation  (32)  states  that  the  deptli  of  the  thermal  layer  should  be  as 
great  as  the  steadv-state  vaUie.  The  second  alternative  of  Eq . (32) 
follov/s  from  Eq.  (23)  and  (26).  By  combining  Eq.  (22)  and  (24)  tlirough 
(32),  Margolin  and  Cbuiko  (1963)  arrived  at  the  following  criteria  for 
ignftion  of  tlic  pore  wall: 

Eroiii  Eq.  (31)  (tt?mperature  criterion) 
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From  Eq.  (32)  (depth  criterion) 


(v  - r)rd 

__S E_ 

4Nua  a 
8 s 


> 1 


(34) 


Both  inequalities  implicitly  contain  the  hydrodynamic  criterion  v 
(see  Eq . (14) ) . 


8 


> r 


It  should  be  noted  that  Eq.  (33)  and  (34)  are  quite  general,  within 
the  limits  imposed  on  their  derivation,  because  Vg  and  Nu  depend  upon 
unspecified  experimental  or  model  conditions.  Margolin  and  Chuiko 
(1965)  presented  Eq.  (33)  and  (34)  in  a schematic  graphical  display  (see 
Figure  50)  of  pressure  vs.  pore  diameter  without  specifying  the  relevant 
experimental  or  model  conditions  and  without  providing  any  details  of 
the  transformations  required  to  obtain  the  graph.  Therefore,  discussion 
of  the  implications  of  Figure  50  must  be  made  in  the  absence  of  any 
resolution  of  the  conditions  of  its  applicability.  Curve  1,  represent- 
ing Eq.  (33),  separates  the  p-dp  plane  into  two  regions,  the  upper  one 
denoting  conditions  for  which  the  surface  is  heated  to  the  ignition 
temperature  or  above.  The  horizontal  nature  of  Curve  1 suggests  that 
pore  diameter  is  not  present,  even  implicitly,  in  Eq.  (33).  Actually  Vg 
and  Nu  may  depend  on  pore  diameter  while  T„  could  be  influenced  both  by 
pressure  and  pore  diameter  owing  to  the  effect  of  pressure  on  flame 
standoff  distance  and  the  consequent  temperature  of  gases  flowing  into 
the  pore.  The  effect  of  the  variables  in  Eq.  (33)  on  the  position  of 
Curve  1 cannot  be  estimated  until  the  dependence  of  (v„-r)/r  on  pressure 
is  established;  this  in  turn,  is  not  a general  result  but  requires 
assignment  of  model  conditions.  Curve  2,  representing  Eq.  (34),  also 
separates  two  regions;  the  lower  consisting  of  combinations  of  pressure 
and  pore  size  for  which  the  thermal  layer  is  too  thin  to  assure  igni- 
tion. The  location  of  Curve  2 is  undetermined  unless  the  pressure  and 
pore  size  dependencies  of  the  terms  in  Eq.  (34)  are  known  or  assumed. 


Curves  1 and  2 delineate  four  regions  in  the  p-dp  plane.  In  Region 
A,  both  conditions  for  pore  wall  ignition  are  satisfied,  resulting  in 
the  onset  of  convective  burning.  In  Region  B,  neither  condition  is 
satisfied,  so  that  combustion  does  not  penetrate  the  pores.  However, 
there  is  some  inflow  of  hot  gases  which  heat  the  pore  walls  above  the 
ambient  temperatures  and  lead  to  slight  augmentation  of  the  conductive 
burning  rate.  Region  C corresponds  to  a hot  enough  but  not  a thick 
enough  thermal  layer.  The  result  is  known  as  forced  pyrolysis,  which 
would  cease  with  removal  of  the  stimulating  source  of  energy.  Region  D 
represents  a condition  in  which  a thick  layer  of  the  pore  wall  is  heated 
to  some  temperature  below  the  ignition  point.  The  result  would  be 
either  a thermal  explosion  (result  of  synergistic  interaction  between 
temperature  and  rate  of  heat  release  by  chemical  reaction)  of  the  ther- 
mal layer  or  a significant  augmentation  of  the  conductive  burning  rate. 
Region  E,  included  in  the  reference,  is  said  to  represent  turbulent 
combustion  of  gases  which  flow  into  the  pores,  but  further  details  are 
sketchy . 
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h LUUKL  5U.  bc'liematic  ot  Poro 
Ignition  Regime's. 


The  foregoing  analysis  and  remarks  are  strictly  appl  ic.abl  (?  to 
smooth  wall  pores.  Brief  mention  of  roughne.ss  of  pore  wails  is  made  by 
Margolin  and  Chuiko  (1965).  A protuberance  may  st'rve  as  an  ignition 
center  if  two  conditions  are  satisfied.  Define 

1 - D^/u  (35) 

p c s 

as  the  characteristic  heating  time  of  the  protuberance  having  character- 
istic dimension  U . The  conditions  are 
c 

r > T > ( (36) 

P a 

where  i and  Tyi  are  given  by  liq.  (24)  and  (27).  The  first  inequality  of 
Eq . (36)  requires  that  ttie  available  heating  time  be  sufficient  to  heat 
the  particle  to  the  ignition  temperature  wlille  tlie  second  inequality 
assures  the  formation  of  a thermal  layer  of  sufficient  thickness,  equal 
to  the  steady-state  value.  If  Eq.  (36)  is  satisfied,  the  ignition 
Region  A of  Figure  50  may  extt'nd  into  Region  D. 

Thun  far,  eonditiou.s  for  ()ore  wall  ignition  by  inflowing  hot  gases 
have  been  described.  Margolin  and  Chuiko  also  present,  without  deri''a- 
t iini , the  conditions  for  ignition  by  other  mechanisms,  Including; 
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1 . 

2. 

3, 

4. 


P n'pagiiL ion  of 


conibus I. ion  into  pores 


filled  with  conibusLible 


gases. 

Propagation  of  combustion  through  flowing  hot  gases. 

Ignition  of  pore  wall.s  by  thormal  ttxplo.sion  of  reactivi'  pore 
gases . 

Ignition  of  pore  walls  by  adiabatic  compression  of  jiore  gases. 


While  these  alternatives  are  viable  meclianism.s  for  pore  wall  igni- 
tion, neither  the  assumptions  needed  for  the  derivations  nor  the  condi- 
tions of  applicability  are  stated  in  tlio  cited  reference. 

Despite  the  generality  of  F.q . (33)  and  (34),  few  applications  havt' 
been  made  to  particular  model  or  experimental  conditions.  Margolin  and 
Chuiko  (1966)  rewrote  F.q.  (34)  as 


where 


4Nu 


V - r d d 

^ -E  > 1 


r 


(37) 


t = u /r  (38) 

s H 

f = a /r  (39) 

g g 

denote  steady-state  thermal  layer  Lbiekaes.s  for  the  solid  and  gas.  For 
the  case  of  .spontaneous  penetration  of  combustion  (pressure  difference 
generated  by  combustion)  the  functional  relationship 


V 

8 


d 


c 


T T T \ 

T.’  T ’ T j 

* o o/ 


(40) 


is  pre.sented,  where  7,  are  other,  unspecified  rbarac  ter  is  t ic  combustion 
zone  dimensions,  Pr^  are  related  Prandtl  numbers,  Le  is  the  hewi.s  num- 
ber, Tj^  is  melting  temperature,  and  the  other  symbols  have  been  previ- 
ously defined.  Assuming 


1.  v^  >>  r 

2.  is  proportional  to  ? 

3.  The  function  F i.s  only  weakly  dependent  on  its  arguiiii-nts 

4.  Convective  burning  is  independent  of  ' 


Margolin  and  Chuiko  (1966)  combine  Eq . (37)  and  (40)  and  arrive  at  Llie 

requirement 


d /C  = e c 1 d /I  = const  =•  An 


(41) 
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E 


I 

I 

I 


Equation  (41),  known  as  the  Andreev  criteric'n,  is  thus  seen  to  be  re- 
lated to  Eq . (34),  which  expresses  the  ignition  requirement  ol  the 
formation  of  a thermal  wave  thickness  of  appropriate  dimensions.  It  is 
noteworthy  that  the  Andreev  criterion  does  not  involve  a surface  temper- 
ature requirement;  hence,  it  would  predict  ignition  of  pore  walls  in 
both  Regions  A and  D (Figure  42).  Tlie  liydrodynamie  criterion  is  satis- 
fied by  assumption  (1);  otlierwise  tlie  Andreev  criterion  has  limited 
applicability. 

In  the  only  other  reference  to  the  general  ignition  conditions, 
Belyaev  (1973),  in  discussing  result.s  of  manometer  bomb  tests  (si’e 
Section  3. 1.2.2  of  this  review),  presents  tlie  following  criterion 
witliout  derivation. 


l+2n  2 

p d = con.stant 

P 


Equation  (42)  can  be  deduced 


from  Eq.  (34) 


by  assuming 


1. 

2. 

3. 

4. 


Vg  >>  r 

Vg  is  proportional 

a„  is  proportional 
A n 


to  r 

to  1/p,  valid  at  low  pressures 


(42) 


Belyaev  (1973)  could  not  account  for  the  break  in  the  slope  of  the 
(critical  pressure)  vs.  d^  curve  witliout  attributing  the  piienomenon  to 
some  failure  of  the  assumptions  made  in  deriving  Eq.  (33)  and  (34). 
However,  compari.son  of  Figures  5 and  30  (sec  Figure  51)  suggests  tliat 
the  break  may  be  explained  by  a cliange  in  the  relevant  criterion.  For 
small  pores  and  iilgh  pressures,  the  experimental  data  fall  on  Curve  2, 
separating  Regions  A and  C (Figure  50);  hence,  the  criterion  is  the 
formation  of  a sufficiently  thick  thermal  layer.  For  large  pores  and 
lower  pressures,  the  data  correspond  to  Curve  1,  separating  Regions  A 
and  D,  where  the  criterion  is  the  attainment  of  the  required  .surface 
Lemperature.  Differeui  break  points  would  be  expected  for  different 
materials  having  differc-nt  properties.  The  data  for  RDX  (Figure  5) 
probably  show  no  break  simply  because  they  do  not  extend  to  large  enough 
pores . 


Godal  (19/0)  jirovides  a qualitative  theory  to  explain  experimental 
re.sults.  By  equating  the  rate  of  heat  evolved  in  rombu.stion  of  the 
walls  of  flat  crack  to  the  rate  of  heat  lost  (to  these  sami'  walls)  he 
arrives  at  the  following  relation: 

t+  = 4A  (T  - T )/re  Q (43) 

* g g s s s 


73 


NWC  TP  6007 


I'lgure  50  (Reprodui-cd)  . 


FIGURE  51.  Comparison  of  Figures  5 and  50. 


where 

= thre.shold  crack  width  for  penetration  of  combustion 
\g  = ga.s  conductivity 
Tg  = gas  temperature 
Tq  =»  propellant  surface  temperature 

r =■  conductive  burning  rate 

pg  =•  propellant  density 
Qg  = heat  of  combustion 

Equation  (43)  exhibits  the  correct  qualitative  relationship  between 
threshold  crack  size  and  burning  rate  and  also  predicts  an  Inverse 

relationship  between  crack  width  and  heat  of  combustion.  However,  it 

also  predicts  an  increasing  threshold  crack  size  with  increasing  gas 
temperature.  An  additional  shortcoming,  stated  by  Godai,  is  the  failure 
to  explain  behavior  of  aluminized  propellants, 

4.1.3  Propagation  of  Ignition  Front  in  Single  Pore 

The  preceding  sections  have  reviewed  the  available  literature 
dealing  with  two  of  the  necessary  coiidition.s  for  transition  from  con- 
ductive to  convective  burning  in  a single  pore;  flow  of  hot  combustion 
gases  into  the  pore,  and  ignition  of  the  pore  walls  by  the  hot  gases. 

The  tliird  requirement  is  that  the  ignition  front  continue  to  propagate 
into  the  pore  at  a rate  exceeding  the  normal  conductive  burning  rate. 

The  only  avnilatile  simplified  model  which  is  related  to  the  phenomenon 
[Krasnov  (1970)]  i.s  depicted  schematically  in  Figure  52.  The  hot  com- 
bustion g.ises  with  t einpe r .a ture  Tj,,  velocity  Vg]),  density  i gp  flow 


I 

/ u 


A 


FIGURE  52.  Simplified  Model  of 
IpniLion  Propagation  into  Single 
Pores . 


into  the  raoutii  of  a eiia'ui.ir  pore  of  diameLet  dp.  After  Ltaversiiig  a 
cooling  section  of  length  gas  flows  away  at  the  lower  tempera- 
ture T.3,  velocity  and  density  Pj-c-  cooling  process  is  assumed 

to  occur  at  a constant  temperature  difference,  as  discussed  in  the 
previous  section  on  ignition.  Using  a coordinate  system  moving  with  the 
ignition  front  with  the  velocity  Vj^g,  energy  and  mass  conservation  for 
the  ignition  section  are  written 


p ,C  (T  - T ) (v  - V , )Trd  /4  = p c (T.  - T )iTd  v,  6 (44) 

^gh  g'  g o-'''  gh  ig^  p ' 'p  p * o"  p ig  t ^ ' 


P ig 


P . (v  - V . ) = D (v  - V . ) 

gh  gh  ig  ge  gc  ig 


(45) 


Equation  (44)  equates  the  heat  lost  by  the  gas  to  the  heat  gained  by  the 
heated  layer  of  propellant,  whose  thickness, 
pared  to  the  pore  diameter.  By  eliminating 
(45) , one  may  obtain 


is  assumed  small  com- 
(44)  and 


V between  Eq 
gi' 


6 

t 


V - V . c ,1  T - T 
8t~  ig  R gr  ^ o 

4v.  e i.  T.  - T 

Ig  s s o 


(46) 


for  the  thickness  of  tlie  tliermal  layer.  The  cooling  length  may  be  found 
f rom 


h(T 

R 


T.)7,L  d 

* c 


I’ 


o C (T 

gr  R g 


T )(v 
o gc 


V . ) Trd 

Jg  P 


/4 


(47) 
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which  equates  the  heat  transferred  throuph  the  gas  film  to  the  heat  lost 
from  the  flowing  gas.  lij/o  conditions  are  considered;  smooth  jinri'  walls 
and  rough  pore  walls  with  protuberances  of  the  optimum  size  u /r  (thick- 
ness of  steady-state  thermal  layer). 

Smooth  Walls 


The  time  available  for  wall  heating  is  llmlt<;d  by  the  time  required 

to  burn  the  distance  L at  the  conductive  rate  r (ef.  Section  4.L.2) 

c 


i = L /r 

t 

A 

The  thermal  layer  thickness  at  time  t is 
^t  = ("s^^ 

Combining  Eq.  (46)  through  (49)  gives  the  result 

no  c A ii’'  - T )^ 

= 1 + « o 


V 

-^'1  = 


V . 

Ig 


Nup  c A (T  - T )(T  - T*) 

g g g g o g 


Rough  Walls 


The  time  required  to  heat  a particle  of  optimum  size  is 
/ 2 

T = a /r 
P 

Combining  Eq , (46)  through  (48)  with  Eq . (51)  gives  the  result 


(48) 


(49) 


(50) 


(51) 


4Nuu  a 

1 + ^ 

2,2 
r d 

P 


(52) 


It  may  be  concluded  that  for  either  smooth  or  rough  wall  pores,  the 
cooled  gas  velocity  exceeds  the  ignition  front  velocity.  For  smooth 
walls,  the  ratio  Vg(,/vj^g  (referred  to  from  here  on  as  ignition  lag) 
should  be  independent  of  diameter,  and  conductive  rate.  IncrL’asing 

the  pressure  should  decrease  the  ratio  since  og  and  Tg  can  only  increase 
with  pressure.  For  rough  walls,  the  ignition  lag  should  he  indcpi'ndent 
of  Vt,;-,  but  should  decrea.se  with  increasing  values  of  pressuri'  and  pore 
diameter. 


■k 

Use  of  Eq . (49)  is  actually  Inconsistent  with  assuming  a constant 
temperature  for  the  thermal  layer  since  Eq.  (49)  is  based  on  a variable 
tempe  rature . 
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Kxpi'r imi'nt.'i  1 rosulLs  (si’e  Soctinn  3.  1.3.4,  ,'itul  Kiguri-  20  siigp,osL 
that  LIilTi’  is  a dopi'iidfiico  nl  ignition  l.ig  on  d^  tor  pores  smnJ  ler  than 
2 mm  which  agrees  qualitatively  with  Kq . (52)  tor  rough  pore  walls.  On 

the  other  hand,  experimental  results  for  larger  pores  appear  to  be  .size 
independent  as  required  by  hq ■ (50)  fot  smooth  pore  walls.  A possible 
explanation  of  this  anomaly  is  that  roughnes.s,  represented  by  a constant 
size  protuberance,  would  be  relatively  more  important  in  a small  channel. 
A seriou.s  unexplained  departure  of  experimental  results  from  thettry  is 
the  dependenee  of  expi'r  iimai  t.i  1 ignitinu  lag  on  v^^. . It  is  thus  concluded 
that  the  claim  of  good  agreement  between  experiment  and  theory  by 
Krasnov  (1970)  is  not  well  supported  by  the  data  presented. 

It  should  be  noted  that  the  theoretical  results  derivi'd  in  this 
section  are  not  depend<.'nt  on  any  particular  model  or  experimental  eondi- 
tions,  such  as  constant  pressure  (prossun'  i.s  not  even  nu'ntioned)  or 
channel  end  condition.  i)n  the  other  hand,  the  experimi-nt  was  carried 
out  in  a very  specialized  const.aat  pressure  apparatus  in  which  the  gas 
flow  into  the  channel  displaced  a liquid  through  a valve  ,jL  the  unigniLed 
end . 


Although  the  paper  [Krasnov  (1970)]  is  relati.'d  to  the  problem  of 
propagation  rate  of  the  ignition  front,  it  tiocs  not  present  an  explicit 
formula  for  v^g,  but  only  for  the  ratio  Vg,./vj^g.  The  requisite  an.alysis 
to  determine  the  cool  gas  velocity  Vg^,  wnich  UL'pends  on  Vg|,  and  thus  on 
model  conditions,  lias  not  been  carritid  out.  Thus,  at  this  timi.',  there 
is  no  available  complete  simplified  analysis  whicli  leads  to  .a  criterion 
for  fulfillment  of  the  third  requirement  for  convective  burning  in  a 
single  pore,  viz.,  v^^  > r. 

4.2  POROUS  BED 

The  simple  theorie.s  of  porou.s  bed  convective  burning  are  directed 
exclusively  at  the  solution  of  the  fluid  dynamic  problem  of  transient 
flow  into  a iiorous  bed  under  several  bonndary  conditions.  Only  the  mas.s 
and  momentum  equations  arc  involved,  while  the  energy  equation  i.s  elim- 
inated by  assuming  isothermal  flow.  This  simplification  alone  precludes 
the  e.stabl  ishment  of  ignition  condition.s  .since  tln-rc  is  no  calcul.ition 
of  pore  wall  lieating  during  the  flt>w.  The  only  attempt  to  consider 
ignition  conditions  is  the  heuristic  extension  of  tin'  Andreev  (.riterion 
to  porou.s  bed.s  by  defining  an  eqniv.ilent  average  hydraulic  pore  diameler 
based  on  permeability  nn-asuremen  t s . The  siun'es.s  of  such  ,i  jiroccdiiri'  has 
already  been  di.s(’ussed  in  Section  3.2.2.  1.  Finally,  there  .ari'  no  sim- 
plified theories  addressing  the  problem  of  propagation  of  the  combiist ion 
front  into  a |)orons  bed. 
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4.2.1  Fluid  Dynamics 

The  fluid  dynamics  associated  with  simplified  modeling  of  convec- 
tive burning  in  porous  beds  is  summarized  by  Belyaev  (1973).  The  analy- 
sis uses  the  momentum  equation  in  the  form  of  Darcy's  law 

u = 1/^^,  dll/dt  = - k/u  3p/3x  (53) 

where 


u = volumetric  flow  rate  per  unit  cross-sectional  area  of  the  bed 
(not  the  pores) 

k = proportionality  constant  called  permeability 
u = fluid  viscosity 
3p/3x  = pressure  gradient 

The  conservation  of  mass  is  expressed  by 


IT  = 3^ 


(54) 


where 


m = porosity 
p = fluid  density 


Tlie  basic  equation,  referred  to  as  the  filtration  equation  by  the 
Russians,  is  obtained  by  combining  Eq.  (53)  and  (54)  to  give 


ilil  = L,  ^ (p  i£) 

3t  um  3x 


(55) 


In  the  derivation  of  Eq . (55)  it  is  assumed  that  k and  u are  constants. 
In  general,  the  energy  equation  and  equation  of  state  would  be  used  to 
eliminate  eitlier  p or  p from  Eq . (55);  however,  in  the  simplified  ap- 
proach, a polytropic  relation 

p = (56) 


is  assumed.  Elimination  of  ()  between  Eq.  (55)  and  (56)  yields  the 
following  equation  for  pressure 


where 


p_m_  J)_l'  _ _L>^P 
nkp  It 

n+1 


(57) 


(58) 
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The  only  known  analytical  solution  of  Eq.  (57)  is  when  p is  nearly 
constant  (small  pressure  gradient),  in  which  case  P satisfies  the  dif- 
fusion equation.  Belyaev  makes  no  further  use  of  this  solution  but 
turns  to  the  solutions  of  Eq.  (57)  for  large  pressure  gradients.  Two 
additional  assumptions  are  made: 


1.  The  initial  pore  pressure  may  be  neglected  in  comparison  to 
the  external  pressure. 

2.  The  process  is  isothermal  (n=l) 


The  ensuing  equation 


ilH 

at 


k 

2um 


(59) 


is  solved  for  an  infinite  porous  medium  and  for  a porous  medium  of 
length  L.  Exact  analytical  solutions  of  the  non-linear  equation  (Eq. 
(59))  are  tiot  known  but  approximate  solutions  may  be  obtained. 


Infinite  Medium.  The  approximate  solution  of  Eq . 
Stant  exLelllal  pressure  p is 

p(x,t)  = 

where 


2.29  C 


2V^Po 


(59)  for  a con- 


(60) 


C 


2 

2 


k/ 2m(i 


(6’  ) 


and  the  pressure  is  seen  to  decrease  linearly  with  increasing  distance. 
The  position  of  the  advancing  gas  front  is  found  by  equating  the  expres- 
sion in  parentheses  tc  zero.  Hence 

X(t)  = 1.62\/"kp^t/mp  (62) 

The  velocity  of  the  leading  edge  of  the  penetrating  gas  is  found  by 
differentiating  Eq . (62)  with  respect  to  time 


v (t)  = 
g 


81  J kp  /mil  t 


(63) 


The  fluid  aynamic  criterion  fot  gas  penetration  (Eq.  (14))  c'an  be  stated 
in  an  alternate  manner:  gas  must  penetrate-  to  a depth  exceeding  the 
steady-state  thermal  layer  thickne.ss  ('a  It)  In  the  time  required  for 
this  laye>r  to  burn  (ujj/r^).  Hence,  substituting  X=r(  /r  and  t=  < /r^ 
into  Eq.  (62)  gives  the  fluid  dynamic  criterion  ' ^ 

kp  /iciu  ' 5 0.  38  (64 ) 

o ,s 
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A numerical  example,  with  k=10  darcy  (10  cm  ),  m=0.05,  ii=3xl0 

poise,  ag=10~3  cm^/s,  and  Pp=50  atm  (5  MPa)  Indicates  that  gases  will, 

in  fact,  penetrate  porous  beds  having  very  low  permeabilities.  Ignition 

is  of  course  not  implied. 

Finite  Medium  of  Length  L.  The  solution  of  Eq . (59)  for  an  imper- 
meable boundary  at  x=L  is  given  by  Eq . (60)  until  the  gas  reaches  the 
boundary.  At  that  instant  (t'),  the  boundary  condition  9p/lx-0  must  be 
imposed,  leading  to  the  approximate  solution 


IPTv  — V 

p(x.i;)  = p^jl — Y — - exp[-  3kp^(t  - t')/muL  ] (65) 

1 ^ I 


At  the  boundary  (x=L),  Eq . (65)  becomes 


p(L,t)  = p^  1 


exp[-  3kp^(t  - t')/mpL  ] 


for  times  greater  than  t'.  The  time  required  for  the  gas  to  reach  the 
boundary  is  found  from  Eq.  (60)  by  setting  x=L  and  p(x,t)=0. 


" ■ 2.62  kp 

An  experiment  was  conducted  to  verify  the  above  analysis.  A propellant 
charge  was  prepared  from  a mixture  of  90%  AP  and  10%  polystyrene. 
Conditions  of  the  test  (pressure=25  atm  (2.5  MPa],  m=0.15,  k=10“3  darcy, 
ij=0.03  cP)  were  such  that  convective  burning  was  precluded.  Measurement 
of  pressure  at  the  closed  end  of  charges  of  lengtlis  1 and  2 cm  confirmed 
that  the  filtration  time  t'  is  given  with  reasonable  accuracy  by  Eq . 
(67).  This  agreement  is  taken  to  justify  the  as.sump'- ions  made  in  deriv- 
ing Eq . (67).  There  are  two  additional  applications  of  the  analysis. 
Both  depe.J  on  the  definition  of  pressure  relaxation  time  for  a closed 
channel  a.s 


In  this  additional  time  after  the  arrival  of  the  gas  at  the  floscd  end, 
the  end  pressure  increa.sus  surii  that 


1 - 1/e  = 0.6: 
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One  application  is  the  statement,  without  further  analysis,  that  the 
presence  of  the  closed  end  restricts  the  inflow  of  gases  and  ma_^  pre- 
vent convective  burning  if  ignition  of  the  pores  Joes  not  occur  before 
pressure  relaxation.  Although  not  st.ited,  this  relaxation  time  would  be 
the  total  time  required  for  the  flow  to  reach  the  end  plus  the  addi- 
tional time  t^.  The  other  application  relates  to  the  application  of  Eq . 
(17)  to  a porous  bed  for  the  case  of  rising  pre.ssure.  Tlie  limitation  of 
Eq . (17)  is  the  equality  of  the  relaxation  time  of  the  bomb  pressure 
(p/p)  and  the  time,  tj-.  Thus,  the  use  of  Kq.  (17)  to  calculate  velocity 
of  the  gases  flowing  into  a blind  porous  lied  is  limited  by 


dt  2 

mpL 


(70) 


-3 

Using  the  values  p=50  atm  (5  MPa),  k=10  darcy,  m=0.15,  L=5  cm,  ii=0.03 
cP  gives  a limiting  value  of  dp/dt  of  50  atm/s. 


4.3  MISCELLANEOUS 

Sections  4.1  and  4.2  have  summarized  the  available  literature' 
dealing  with  theoretical  approaches  to  early  stages  of  t.-onvective  burn- 
ing in  single  pores  and  porous  beds.  This  final  section  presents  .sev- 
eral tlieoretical  results  which  do  not  fit  the  fluid  dynamic,  ignition, 
or  propagation  requirements  of  convective  burning.  Nevertheless  they 
are  a part  of  the  body  of  theoretical  knowledge  and  are  included  for 
completeness . 

4.3.1  Effect  of  Mel t Layer 

The  issue  of  cf fectivene;  s of  the  melt  layer  in  hindering  the 
development  of  convective  burning  is  still  controversial.  IVo  theoreti- 
cal efforts  to  approach  the  problem  have  been  reported  by  l^elyaev  (1973). 
In  the  first,  the  solution  of  the  s teady-.state  heat  conduction  equation 
for  a moving  solid  is  applied  to  tlie  calculation  ol  the  melt  layer 
thickness 


^me  i t 


A 

I'cr 


m o 


whe  re 


„meLt 


Tx 


= projierties  of  the  pielt 
= melting  temperature 
= initial  propell.int  temperature 
= critical  temperature 


(71) 
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Although  not  defined,  Is  probably  the  highest  temperature  at  whleh 
the  melt  can  exist  either  because  of  boiling  or  di'composition.  Melt 
layer  thicknesses  v/ere  calculated  from  E (71)  (.parameter  values  not 
provided)  and  compared  with  experinentell^  determined  critical  pressures 
for  breakdown  of  conductive  burning.  The  experimental  conditions  in- 
volved blind,  non-embedded  porous  beds  with  diameters  of  10  mm,  lengths 
ranging  from  40-70  mm,  permeability  of  10“^  darcy,  and  loading  density 
of  0.05  g/cm^  (HLDCB) . Results,  shown  in  Table  8 indicate  that  break- 
down pressure  increases  with  calculated  melt  layer  thickness.  Moreover, 
the  noil-fusible  materials  such  as  NC,  AP-based  composite  propellants  (AP 
Is  regarded  as  non-fusible  even  though  contrary  evidence  exists 
[Hightower  and  Price  (1967  and  1968)  and  Boggs  and  Kraeutle  (1969)], 
mercury  fulminate,  and  lead  azide,  exhibit  lower  critical  pressures  than 
shown  in  Table  8.  It  is  concluded  that  the  melt  layer  Is  effective  in 
impeding  convective  burning. 


TABLE  8.  Comparison  of  Breakdown  Pressure  and 
Melt  Layer  Thickness  for  Fusible  Substances. 


Melt  layer  tbi  ckne.s.<? , pm 


Subs  tauce 

100  atm 

10  (MPa) 

300  atm 

30  (MPa) 

Critical 

atn 

pressu 

(MPa) 

TNT 

50 

18 

2,000 

(200) 

Picric  acid 

35 

1 

i2 

800 

(80) 

PETN 

13 

3 

550 

(55) 

RDX 

5 

2 

2 50 

(25) 

The  second  theoretical  approach  involving  the  melt  layer  states 
without  derivation  that  the  thre.shold  burning  rate  of  a fu.sihle  pm- 
pellant  is  given  by 

12(1  - m) (p  - p' )r 


whete 


e = solid  density 
p ' = melt  density 
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= conducti\'e  burning  rate  at  a prossura  of  1 atm  (0.1  MPa) 
m = porosity 

R = characteristic  particle  size  (radius?) 

No  application  is  made  to  evaluation  of  experimental  data;  indeed,  it  is 
not  even  clear  which  type  of  experiment  is  relevant  to  Eq . (72).  A 
final  comment  is  made  that  in  evaluating  the  Andreev  number  for  fusible 
propellants,  the  pertinent  dimension  to  use  is  the  melt  layer  thickness. 

4.3.2  Dynamic  (Combustion  Generated)  Pressure 

In  the  absence  of  a rising  bomb  pressure  or  of  an  iiiqjosed  initial 
pressure  difference  between  bomb  and  pore,  penetration  of  combustion  is 
dependent  upon  the  usually  small  pressure  difference  generated  by  the 
combustion  process.  This  pressure  gradient  is  needed  to  cause  the  flow 
of  gases  away  from  the  burning  surface  to  the  surrounding  atmosphere. 

The  derivation  (Belyaev  (1973)  and  Bakhman  (1965)]  is  based  on  conserva- 
tion of  mass  and  momentum. 

Mass ; 


p r = p'v'  - p,v, 

P 1 1- 

Momentum; 

p'  + (I  rv'  = p,  + p rv, 

P 1 P 1 

where  the  primes  refer  to  conditions 
the  subscript  1 refers  to  conditions 
Tlie  solution  of  Eq . (73)  and  (7<!)  fo 


(73) 


(74) 

in  tiie  gas  phase  at  the  surface  and 
of  the  final  combustion  products. 
APj  = P'  - Pj^  gives 


Ap,  = (p  r)^  (1/p.  - 1/p')  (75) 

d p 1 

which  becomes,  assuming  Apj<<Pj^  and  using  tlie  perfect  gas  equation  of 
s tate : 


(76) 


where  M is  molecular  weight.  Tlie  dynamic  pressure  increa.scs  with  the 

square  of  the  mass  burning  rate,  but  even  so,  amounts  to  only  a few  torr 

for  most  combustibles.  Mercury  fulminate,  with  a burning  rate  (ppt)  of 

5.9  g/cm3-s,  gener.'ites  a Apj  (jf  13  torr  while  lead  styphnate,  with 

burning  rate  tjf  lOO  g/rra^-s  detonates  upon  ignition.  Little  use  is  made 

of  Eq,  (75)  since  the  pre.ssure  differentials,  even  in  LLDCBs,  exceed  the 

usual  values  of  Ap , . 

d 
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4.3.3  Effect  of  Pressure  Oscillations 


The  possiblity  that  combustion  oscillations  (as  from  combustion 
instability)  could  aid  the  onset  of  porous  burning  in  a porous  bed  was 
investigated  theoretically  by  Margolin  (1961).  The  analysis  is  based  on 
an  assumed  polytropic  flow  using  the  equation  of  mass  and  momentum 
conservation . 


Momentum: 


P 


9v 

at 


+ pv 


3v 

3x 


3p  mtiv 

3x  k 


(77) 


Mass : 


It  + It  <»'')  ■ 0 

where  v is  gas  velocity  and  the  other  symbols  have  their  usual  meanings. 
For  boundary  conditions  of  p=Po  at  x=-°’  and  p=Po+Pl  cosojt  at  the  burning 
face  of  the  porous  bed,  the  results  may  be  summarized  as  follows: 

Conditions  for  onset  of  flov; 


(79) 


2 2 2 
where  N=2u)a  k/r  mp  for  isentropic  flow  and  N=Pq  koi/r  mp  for  iso- 
thermal flow.  In  Eq.  (79),  a is  sonic  velocity,  m is  frequency,  and 
is  average  gas  density-  Conditions  for  penetration  to  the  depth  S=2ti(v- 
4)  /u) 


(80) 


Minimum  pressure  amplitude  and  optimum  frequency  (w^^)  for  penetration  to 
(a)  depth  S 


0)^  = 2iir/S 


^o/mln 


(81) 


(b)  depth  equal  to  thermal  layer  thickness  a /r 
0 2 . 

iji,  = 2 nr  /tt 
i P 


(82) 
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jl  mi; 


2 up  k. 


(83) 


In  Eq . (81),  N(ijjj)  is  the  appropriate  function  N tisentropic  or  iso- 
thermal) evaluated  at  the  frequency  co^. 

No  experimental  verification  of  the  theory  has  been  found;  liowever, 
owing  to  the  occurrence  of  combustion  instability  in  many  operational 
rocket  motors,  it  would  be  interesting  to  investigate  the  possible 
relation  to  DDT  through  assisted  convective  burning. 


5.0  CLOSURE 


The  previous  sections  summarize  the  literature  describing  the  onset 
and  development  of  convective  combustion.  This  final  .section  is  devoted 
to  general  comments  on  the  scope  and  character  of  the  surveyed  liter- 
ature, an  enumeration  of  aspects  of  propellant  systems  that  favor  devel- 
opment of  convective  burning,  and  identification  of  problem  areas  demand- 
ing further  study. 


5.1  OVERVIEW  OF  THE  SCOPE  AND  CHARACTER  OF  SURVEYED  LITERATURE 

Soviet  investigators  have  worked  for  several  decades  on  the  prob- 
lems associated  with  early  stages  of  transition  from  conductive  to 
convective  burning.  Their  investigations  deal  largely  with  model  samples 
having  characterizable  defects  rather  than  actual  propellants  with  more 
random  defects.  Their  sample  Ingredients  and  characteristics  were 
chosen  primarily  to  facilitate  studying  basic  principles  and  not  nece.s- 
sarily  to  optimize  the  mechanical  or  ballistic  properties  of  operational 
rocket  propellants.  As  explained  in  Section  2.1,  the  defect.s  are  char- 
acterizable as  single  pores  with  definite  geometrical  dimensions  or  as 
porous  beds  with  specified  porosity,  permeability,  pore  size,  particle 
size,  .and  other  statistical  quantities.  Porosities  encountered  are 
often  Higher  than  those  that  might  be  encountered  with  propellants. 
Polymeric  materials  appe.ar  to  be  used  only  as  fuels  with  no  attention 
given  to  their  binding  properties  (fuel  beads  mixed  with  oxidizer  parti- 
cles) . 

Experiments  using  the  sample/defect  combination  just  described  lead 
to  results  of  a fundamental  nature,  but  no  means  are  provided  for  ex- 
tending these  results  to  operational  propellants  and  operational  condi- 
tions. Within  the  Limitations  imposed  by  the  sample/defect  combination, 
the  scope  of  the  Soviet  studies  is  rather  broad,  covering  many  combina- 
tions of  co'idition.s . However,  this  coverage  is  not  coni|5iete,  as  is 
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evident  from  an  examination  of  the  moriihologica  1 chart  (Table  1,  I'.age 
9);  for  example,  littlf*  work  has  been  reported  for  the  open/open  samplL- 
in  the  HLDCB . In  addition,  there  is  an  apparent  lack  of  systematic 
approach.  ExampLes  of  this  are: 

1.  There  is  little  effort  to  correlate  experimental  results  with 
tlieory.  Qualitative  agreement  is  sometimes  c-.laimed,  but  no 
examples  were  found  in  which  experimental  and  tlieoret  i cal 
findings  were  compared  graphically  or  tabularly. 

2.  There  is  no  presentation  of  experimental  data  in  which  a given 
model  sample  is  subiected  to  a spectrum  of  experimental  condi- 
tions. Claims  that  this  was  done  are  supported  cntly  by  quali- 
tative statements  as  to  tlie  similarit  it's  or  differences  of 
resul ts . 

3.  In  the  case  of  different  model  samples  under  the  same  condi- 
tion.s,  the  situation  is  better,  but  it  is  somctimi.-s  difficult 
to  determine  the  conditions.  An  outstanding  example  is  the 
frequent  failure  to  report  whether  or  not  a length  effect 
exists  for  a blind  porous  bed  or  single  pore.  Since  many  of 
the  test  samples  are  at  least  geometrically,  If  not  fluid- 
dynamicaliy  blind  (constructed  by  pouring  or  pressing  granular 
ingredients  into  a dead-end  container),  this  is  a serious 

def iciency . 

Most  of  the  studies  are  addressed  to  the  determination  of  the 
conditions  required  for  the  onset  of  convective  burning,  where  the  onset 
is  identified  by  the  appearance  of  a singularity  in  the  experimental 
results,  e.g.,  a break  in  the  pressure  vs.  time  curve.  More  elaborate 
experiments,  designed  for  determining  the  rate  of  propagation  of  the 
convective  combustion  front,  have  been  reported  only  for  the  Ll.DCB . 
Results  of  both  types  of  experiments  are  generally  shown  in  graphical 
form  with  an  occasional  empirical  correlating  equation. 

A final  comment  concerns  Che  manner  of  presentation  in  the  Russian 
papers.  While  the  language  barrier  may  be  involved,  the  style  seems 
unusually  terse,  often  sketchy,  wltli  insufficient  information  to  su[)porL 
the  conclusions  presented  or  to  enable  Che  reader  to  draw  the  same 
conclusions.  In  certain  instances,  ail  relevant  lactors  have  not  been 
considered,  leading  to  unre.solved  controversies.  An  example  is  the 
di-screpancies  found  for  the  effect  of  the  melt  layer  on  inhibiting 
convective  burning.  Here  the  effect  of  rate  of  change  of  pressure  is 
completely  ignored  and  no  consideration  is  given  to  melt  proiierties 
other  than  melting  point.  A further  difficulty  is  the  lack  of  documen- 
tation of  controversy.  Many  authors  not  only  fall  Co  reference  the 
work  cf  otihers  that  conflicts  with  tiielr  own,  but  in  several  cases  an 
individual  has  presented  e:vldence  that  directly  contradicts  results  he 
previously  presented,  yet  no  mention  ts  made  of  the  disparity. 
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5.2  ATTRIBUTES  FAVORING  CONVECTIVE  BURNING 


As  described  earlier  (Section  2.2),  there  are  general  requirements 
for  the  breakdown  of  conductive  burning  of  a sample.  These  are,  assum- 
ing the  existence  of  a defect(s). 


1.  Flow  of  gases  into  the  defect  or  porous  bed  I (onset  of  loii- 

2.  Ignition  of  pore  walls  j vective  burning) 

3.  Propagation  of  the  ignition  front  at  a rate  greater  than  the 
conductive  rate  (acceleration  of  convective  burning) 


Tlie  above  sequence  describes  ignition  occurring  as  a result  of  flow  of 
hot  gases  (ignores  effects  of  adibiatic  compression,  thermal  explosion 
of  pore  gases,  and  propagation  of  combustion  into  reactive  pore  gases) 
As  one  means  of  summarizing  the  literature  survey,  we  list,  with  bric;f 
comments,  those  attributes  of  propellants  which  provide  potential  for 
convective  burning  through  tlielr  relation  to  one  or  more  of  the  above 
general  requirements. 


Attributes  Favoring  Gas  Flow  Into  Defect 


1.  The  open-open  configuration  provides  a flow  patli  for  the 
initial  gases  to  be  displaced  by  the  inflowing  hot  gases. 

2.  Gas  may  flow  more  readily  into  a flat,  blind  crack  than  into  a 
round  blind  hole  because  of  recirculation  which  provides  an  exit  for 
initial  gases. 


3.  Flow  into  the  defect  is  facilitated  by  the  embedded  sample 
because  of  the  difference  between  bomb  and  defect  pressures  at  the 
Instant  the  conductive  burning  front  encounters  the  defect. 


4.  High  rates  of  pressuie  rise  in  the  bomb  lead  to  pressure 
differentials  which  cause  gases  to  flow  into  defects. 

5.  There  is  conflicting  evidence  concerning  the  effectiveness  of 
a surface  melt  layer  in  preventing  combustion  gases  from  flowing  into 
defects.  Available  data  suggests  that  the  melt  layer  is  more  effective 
under  nearly  constant  pressure  conditions  (LLDCB)  than  under  c:ondltionn 
of  rising  pressure  (HLDCB) . In  regard  to  composite  formulations,  a con- 
troversy exists  as  to  whether  oxidizer  or  fuel  melting  is  more  Important 

6.  Large  diameter  particles,  large  diameter  pores,  and  high 
porosities  lead  to  high  permeabilities  to  gas  flow;  all  favor  the  inflow 
of  gases  but  none  ]jrovldes  a unique  c.urrelation  for  detenu  i nal  ion  ci! 
onset  of  convcc:tlve  burning. 
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A t tribuUtis  Favoring’  I{;niLloii  ol  Puru  Walls 

1.  The  pore  walls  of  an  embedded  sample  are  more  readily  ignited 
because  the  initial  cold  pore  j^ases  are  at  low  pressure  and  provide  less 
dilution  of  the  inflowing  hot  gases  than  would  be  the  case  lor  an  open 
sample  tested  at  high  pressure. 

2.  Pore  walls  are  more  readily  ignited  by  higher  temperature 
inflowing  gases.  Since  these  inflowing  gases  originate  close  to  the 
propellant  surface,  it  is  essentially  the  temperature  profile  in  the  g.is 
phase  which  determines  the  temperature  cf  the  igniting,  hot  gases. 

Titus,  it  is  important  to  consider  the  iinal  flame  temperature  and  the 
flame  standoff  distances  (determined  by  gas  phase  reactions  kinetics). 
Both  features  favor  development  of  convective  burning  in  composite 
propellants  as  compared  to  double-base  propellants.  High  flame  temper- 
ature alone  is  not  sufficient,  as  is  shown  by  the  rsulLs  of  tests  with 
AP/sucrose  mixtures.  lu  a series  of  experiments  using  different  propor- 
tions of  ingredients  to  provide  widely  varying  flame  temperatures,  no 
difference  was  observed  in  the  propensity  toward  convective  burning. 

i.  Materials  of  high  conductive  burning  rate'  gener.aily  tr.ansit 
more  readily  to  convective  burning,  other  conditions  being  the  same. 

This,  in  part,  is  reflected  by  the  Andreev  criterion,  which  represents 
the  requirement  for  establishing  a thermal  wave  of  the  required  thick- 
ness in  the  pore  wall.  Faso  of  propellant  igni tability , as  determined 
in  arc- image  or  laser  ignition  tests,  is  also  expected  to  favor  flash- 
down  into  defects.  Experimental  results  [Derr  and  Fleming  (1973)]  shov; 
a high  degree  of  correlation  between  high  burning  rate  and  ease'  of 
igni tability . Thus  tliere  is  a strong  suggestion  that  similar  factors 
govern  burning  rate,  ignition,  and  fiasbdown.  The  Andreev  criterion 
also  predicts  that  ignition  should  be  favored  by  low  thermal  diffusivi- 
ties  (Low  thermal  conductivitie.s  and  high  heat  capacities)  but  systematic 
experimental  results  are  lacking. 

4.  Rough  pore  walls  may  ignite  more  readily  than  smooth  ones 
because  protrusions  can  serve  as  foci  for  hot  spot  development.  Tiieo- 
reticaily  there  is  an  optimum  roughness  dimension  correspond ing.  to 
thermal  wave  thickness.  Smaller  particles  may  experience  a Lemperalure 
rise  but  l.aek  the  thermal  eapaeity  to  igrite  the  propellant  substrate. 
Larger  particles  fall  to  heat  sufficiently  because  of  their  higher  heat 
capacity  and  higher  heat  losses  to  Llic  adjoining  propellant. 

Att  ributes  Favoring  Accelerati on  of  the  Convective  Front 

Since  the  convective  combustion  front  can  only  evolve  following 
fulfillment  of  the  first  two  requirements,  all  tlie  conditions  just 
enumerated  favor  acceleration.  Additional  observations  arc; 
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L.  Two  opposing  burning  propoiiant  surt'aoos,  as  con  t ras  Li-’d  to  bo 
single  propellant  surface  opposite  an  inert  surface,  contribute  to 
acceleration  of  convective  burning.  On  the  other  hand,  it  is  interesting 
to  note  that  onset  of  convective  burning  is  the  same  Cor  the  two  geom- 
etries. 

2.  The  acceleration  phase  of  convective  burning  occurs  more 
readily  in  longer  channels. 

3.  A pressure  difference  between  bomb  and  defect  may  result  in 
acceleration  of  convective  burning.  As  noted  by  Taylor  (i962a),  very 
small  pressure  differentials  can  result  in  appreciable  changc-.s  in  the 
rate  of  advance  of  the  combustion  front. 


5.3  GENERAL  COMMENTS 

It  is  important  that  many  of  the  requirement.^  for  breakdown  of 
conductive  burning  are  to  be  found  in  the  high  energy  composite  propel- 
lants currently  in  use.  The  heterogeneous  nature  and  high  solids  load- 
ing of  these  f uel/oxidlzer  mixtures  means  that,  unless  perfet;L  bonding 
is  attained  between  the  components,  there  are  incipient  flaws  already 
built  in  which  can  lead  to  convective  burning  under  appropriate  condi- 
tions. Further,  the  gas  phase  reactions  of  these  propellants  are  com- 
pleted at  much  higher  temperatures  and  closer  to  the  surface  than  is  the 
case  for  homogeneous  propellants  (the  traditional  double-base  propel- 
lants). Indeed,  comparative  data  show  the  greater  susceptibility  of 
model  composite  propellants  to  convective  burning  at  pressures  and 
porosities  likely  to  be  encountered  in  operational  situations. 

Studies  on  burning  of  strained  composite  propellants  show  important 
differences  in  burning  mode  when  the  incipient  defects  are  opened  to  the 
flow  of  hot  gases.  Tliis  suggests  tliat  care  be  exercised  in  interpreting 
the  results  of  tests  on  damaged  propellant  if  conditions  leading  to  the 
damage  are  different  from  conditions  of  the  test.  Otherwise,  closed 
bomb  tests  may  fail  to  reveal  a tendency  toward  convective  burning,  not 
because  of  absence  of  defects,  but  because  of  absence  of  adequate  flow 
channels  under  closed  bomb  condltion.s. 


5.4  PROBLEM  AREAS 

Then:  arc  several  general  problem  areas  connected  witli  early  stages 
of  convective  burning.  As  previously  t'mphasized,  availabli’  Soviet  work 
in  the  field  has  dealt  mainly  with  model  systems  (pressed  mixtures  ot 
granular  .solids).  Although  basic  prini:iplen  have  evolved  from  these 
studies  (see  Section  5.2),  tin-  following  ajipcar  to  bo  important  topics 
for  further  study; 
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1.  Quantitative  effect  of  rate  of  pressure  rise  on  onset  of 
convective  burning 

2.  Resolution  of  the  question  of  the  effect  of  a melt  layer. 
Properties  other  than  melting  point  (c.g.,  boiling  point, 
reactivity,  product  accumulation,  thickness  of  melt,  viscosity 
and  surface  tension  of  melt)  should  be  considered. 

3.  Establishment  of  mixture  laws;  i.e.,  how  are  onset  conditions 
for  mixtures  related  to  onset  conditions  for  individual  com- 
ponents . 

4.  Effects  of  defect  geometry,  particularly  length  and  end  condi- 
tion (open  or  blind). 

A second  set  of  problems  currently  under  investigation  at  NWC  is 
related  to  real  propellants  with  both  characterixable  and  random  defects. 
Characterizable  defects  (open  and  blind  round  holes)  have  been  the 
subject  of  a systematic  series  of  studies  [Prentice  (1962,  1977)]  to 
determine  conditions  for  flame  prop.agation.  Combustion  in  r.andom  defects 
is  the  subject  of  an  investigation  by  Boggs  (1976)  in  the  strained 
propellant  studies. 

A third  important  area  concerns  questions  ot  a more  basic  nature 
such  as: 

1.  Methods  of  flame  zone  modification.  In  particular,  are  ad- 
ditives available  which  will  delay,  witb.out  preventing,  final 
gas  phase  reactions  so  that  cooler  gases  will  flow  into  avail- 
able defects? 

2.  Ignition  properties  of  propellants  including  the  measurement 
of  thermal  conductivity  and  specific  lieat  of  components  and 
mixtures . 

3.  Determination  of  interactions  among  ingredients,  i.e.,  the 
lowering  of  decomposition  temperature  of  HMX  by  admixture  with 
AP.  Particle  size  and  pressure  effects  should  be  determined 
on  this  and  other  candidate  ingredients. 

A major  problem  facing  the  investigation  cf  propellants  containing 
random  deiecis  is  the  absence  of  a means  of  uniquely  characterizing 
these  defects.  Equally  challenging  is  the  need  to  relate  such  a char- 
acterization, obtained  under  static  conditions,  to  the  nature  of  the 
defect  existing  under  actual  dynamic  conditions  in  an  operational 
rocket  motor. 
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NOMENCLATURE 


tant  in  Eq.  (11) 
i-sectlonal  area  of  porous  bed 
number  * pCgdr/X 

etrlc  parameter:  0 for  rectangular  slits,  1 for  round  holes 
fraction  of  particles  with  diameter  D 
tant  In  Eq.  (11)  ^ 

Lflc  gas  constant  ■ R/M 
tant  In  Eq.  (13) 
tant  In  Eq.  (56) 
tant  In  Eq.  (61) 

Lflc  heat 
tant  In  Eq.  (13) 

acterlstlc  dimension  of  surface  roughness 
icle  diameter 
- diameter 

hydraulic  diameter 

r hydraulic  pore  diameter  based  on  mass  averaged  particle 
aeter  (Eq.  (8)) 

I hydraulic  pore  diameter  based  on  harmonic  mass  averaged 
llcle  diameter  (Eq.  (9)) 

laullc  pore  diameter  based  on  permeability  measurement  (Eq.  (10)) 

I coefficient  of  heat  transfer 

leablllty  of  porous  bed  to  gas  flow 

Ith  of  propellant  pore  or  porous  bed 

|th  of  cooling  section  in  pore  flow 

Es  number 

Ikness  of  heated  propellant  layer  beneath  pore  wall  at  time,  t 
|idy  -state  thickness  of  leated  gas  adjacent  to  propellant  surface 
^acterlstlc  dimension  of  1th  combustion  zone 

idy-state  thickness  of  heated  propellant  layer  beneath  pore  wall 
Icular  weight 
aslty  =1-6 
I burning  rate 

S flow  rate  of  gas  into  pore 
I accumulation  rate  % gas  In  pore 

L. . J J 1 _ — £ ■ ^ ^ T 


ameter  used  in  analysis  of  oscillating  pressure  [Eq.  (79)  et  seq] 
^elt  number  | 

isure  exponent  in  conauctive  rate  "law".  Poly tropic  exponent  In 


ssure  at  closed  end  of  porous  bed 
tlal  pressure  external  to  porous  bed 
mum  pressure  developed  in  a propellant  defect  during  combustion 
c pressure  rise 
t of  combustion 

e of  heat  loss  from  pore  gas  to  pore  wall 
e of  heat  transfer  across  gas  film  in  pore 


A^rnilable  C> 


0^ 


R 

Re 

Rv 

r 

T 

Tm 

To 

Tw 

T* 

t 

t» 

tr 

U 

u 

V 

^ig 

w 

X 

X 


Universal  gas  constant 
Reynolds  number 

X' 

Conductive  burning  rate,  surface  regression  rate 

Temperature 

Melting  temperature 

Initial  solid  temperature 

Pore  wall  temperature 

Pore  wall  Ignition  temperature 

Time 

Time  required  for  gas  flow  to  penetrate  to  closed  end  of  porous  bed 

Pressure  relaxation  time  of  porous  bed 

Total  voliime  of  fluid  flowing  through  a porous  bed 

Volumetric  flux  density  In  flow  through  porous  bed 

Linear  gas  velocity 

Linear  propagation  rate  of  Ignition  front  Into  pore  or  porous  bed 

Smallest  dimension  between  opposing  propellant  surface  or  between 

propellant  and  Inert  surfaces 

Position  of  gas  front  In  propellant  defect 

Distance  from  entrance  of  pore  or  porous  bed 


Symbols 

a Thermal  diffuslvlty  * X/pc 

3 Define  by  £q.  (30) 

(5  Relative  density  ■ p/Ptmd 

5^  Thickness  of  solid  thermal  layer  assumed  to  be  at  constant  tempera- 
ture (not  the  same  as  ilg) 

X (^efficient  of  thermal  conductivity 

V Dynamic  viscosity  coefficient 

p Density 

Ptmd  Theoretical  maximum  density;  density  of  substance  having  zero  porosity 
T Time  Interval  available  for  hot  gases  to  heat  pore  wall 
Tp  Characteristic  heating  time  of  particle 

Ts  Time  required  to  reach  steady-state  temperature  distribution  in 
pore  wall 

Time  required  to  heat  pore  wall  to  Ignition  temperature 
(j>  “ Ps^*^p  single  round  pore 

» Ps^dij  for  porous  bed  where  d,  is  hydraulic  diameter  (see  d. , 
d2,  and  d3) 

u Angular  frequency  of  oscillating  pressure 

Subscripts  (except  where  otherwise  noted) 
c Cooled  gas 

g Gas  phase 

h Hot  gas 

o Initial 

p Particle 

s Solid  phase 

* Threshold  value 
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